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Abstract 

Correlations play a crucial role in the nuclear many-body problem. We give an overview of 
recent developments in nuclear structure theory aiming at the description of these interaction- 
induced correlations by unitary transformations. We focus on the Unitary Correlation Operator 
Method (UCOM), which offers a very intuitive, universal and robust approach for the treatment of 
short-range correlations. We discuss the UCOM formalism in detail and highlight the connections 
to other methods for the description of short-range correlations and the construction of effective 
interactions. In particular, we juxtapose UCOM with the Similarity Renormalization Group (SRG) 
approach, which implements the unitary transformation of the Hamiltonian through a very flexi- 
ble flow-equation formulation. The UCOM- and SRG-transformed interactions are compared on 
the level of matrix elements and in many-body calculations within the no-core shell model and 
with Hartree-Fock plus perturbation theory for a variety of nuclei and observables. These calcu- 
lations provide a detailed picture of the similarities and differences as well as the advantages and 
limitations of unitary transformation methods. 
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1. Introduction 

Recent years have seen substantial progress in theoretical methods describing the many-body 
problem of low energy nuclear structure in an ab initio sense. Ab initio means from the beginning, 
without further assumptions or uncontrolled approximations. At present for a system of particles 
interacting by the strong interaction, the most elementary degrees of freedom are considered to 
be the point-like quarks and gluons, whose dynamics is derived from the Lagrangian of Quantum 
Chromo Dynamics (QCD). However, in the low-energy regime QCD cannot be treated by per- 
turbation theory because of the confinement phenomenon. The lowest bound states of QCD are 
baryons and mesons, which are the natural degrees of freedom at low energies. In a system at 
very low energy, below about 50 MeV per baryon, only the proton and neutron are left. Other 
baryons like the A-resonance appear only as intermediate virtual excitations. Also the mesons, 
among which the pions are the lightest ones, do not occur as real particles on their mass shell, but 
may be regarded as bosons mediating the interaction among the baryons. 

Recent QCD lattice simulations PJ give hope that in the near future the baryon-baryon inter- 
action might be "measured" on the lattice with sufficient precision. This would complement the 
scattering data, which can give only indirect information on the interaction in form of the phase 
shifts measured at large distances, where the particles do not interact anymore. 

For proton and neutron precise scattering data exist that allow to fit the parameters of nucleon- 
nucleon potential models. Using global symmetries the form of the potential can be written as a 
sum of central, spin-orbit, and tensor interactions. Precise fits to scattering phase shifts reveal that 
the nucleon-nucleon potential cannot be assumed to be local in coordinate space, which means the 
potential depends not only on the relative distance and the spin orientations but also on the relative 
momentum. The minimal momentum dependence besides the spin-orbit terms is of type (as 
in the Argonne VI 8 potential [Q]) or terms like 4^V(f) + V(r)^^ (as in the Bonn A/B potentials 
Recent developments use the approximate chiral symmetry in the light quark sector of QCD 
to establish for the nucleon-nucleon interaction a perturbation scheme in terms of diagrams con- 
tributing to the so called chiral potentials iflEiH]. By adding higher-order terms to the expansion 
they can be improved systematically. However, with increasing number of contributions more 
and more parameters have to be fixed by experimental data. One great advantage of the chiral 
potentials is that two- and three-body forces are treated on equal footing. 

The nuclear structure problem consists in solving the stationary many -body Schrodinger equa- 
tion for A nucleons 

H = E„ \^„) (1) 

and in investigating the properties of the A-body eigenstates by computing observable quan- 
tities as expectation values or transition matrix elements. The subset of the discrete eigenvalues 
E„ represents the excitation spectrum of a nucleus and can be compared to data directly. The con- 
tinuous part of the spectrum corresponds to scattering states which usually have a rich resonance 
structure. 

As the Schrodinger equations already indicates, solving the many -body problem requires two 
main ingredients, the Hamiltonian H and a representation for the many-body eigenstates |^„) or a 
basis in the A-body Hilbert space. 
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For more than about 6 particles Slater determinants |0[y]) provide the most convenient basis 
for treating the many-body problem numerically. They are usually constructed as eigenstates of 
a one-body mean-field Hamiltonian Hq. Matrix elements of one-, two-, or three-body operators 
in A-body space can be calculated easily using an occupation number representation. Modem 
computers and numerical techniques can tackle A-body Hilbert spaces dimensions up to about 
10^*^. These numbers are impressive but by far not sufficient. The dilemma is that realistic nuclear 
potentials, which fit the scattering data at low energy and induce the high momentum components 
observed in bound states of nuclei, possess a strong repulsion at short distances and a strong 
tensor force. Figure |2] illustrates that for the Argonne VI 8 potential. These properties of the 
nuclear Hamiltonian induce short- and long-range correlations in the many-body eigenstate 
Especially the short-range repulsive and tensor correlations cannot be properly represented by 
Slater determinants. The reason is that a Slater determinant is an anti symmetrized product of 
single-particle states antisymmetrization operator): 

, ^2, • • • , ^aIOm) = ^ 0v,(^l) (pvUl) . . . (/>vMa), = (Xi, CTi, Td (2) 

describing the motion of independent particles and as such cannot describe correlations in the 
relative distances (Xi - xj) between particles. Of course. Slater determinants form a complete 
basis and thus can in principle represent any state. However, as we show in Sec. |4]the number of 
Slater determinants needed to describe the ground state that includes the short-range correlations 
exceeds soon any numerically tractable number for A > 4. The Hamiltonian represented in Slater 
determinants results in matrices with large matrix elements even far off the diagonal. In physical 
terms, in this basis the Hamiltonian scatters to very high lying Hq eigenstates. Therefore, one 
introduces an effective Hamiltonian based on the following concepts. 

Common to many conventional derivations of effective Hamiltonians is that one first divides 
the many-body Hilbert space into a so called P-space or model space and a Q-space. The projection 
operators P and Q = 1 - P are defined with the eigenstates of a Hamiltonian Hq, which is usually a 
one-body operator like kinetic energy or the harmonic oscillator Hamiltonian. One is then looking 
for an effective Hamiltonian Heff that should have the same eigenvalues as the full Hamiltonian 
H 

Heff |(D„) = E„ |0„) (3) 

and its eigenstates should be contained in the model-space, i.e. they can be represented as a finite 
sum of Ho eigenstates that span the P-space, 

P-space 

|0„)=P|0„)= l'^[v])(1>[v]|1>.) ■ (4) 

M 

Therefore, Hgff should not connect P- and Q-space, PHgffQ = 0. Usually one is satisfied with a few 
low lying eigenvalues and eigenstates of Eq. @ and regards the higher eigenstates as orthogonal 
rest that carries less and less physics. Therefore, in practice only the decoupling of the lowest 
eigenstates is needed, i.e. (0„|HeffQ = for few low lying £„. A sharp separation between P- 
and Q-space, as given in the mathematical definition of Eq. (H)), is difficult to identify on physical 
arguments anyhow. 
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As the P-space consists of many-body product states that cannot describe short-range correla- 
tions or high relative momenta q, it is regarded as a low-momentum Hilbert space. Sometimes one 
also speaks of a low-energy Hilbert space, but here energy refers to Ho and its eigenvalues and not 
to the true Hamiltonian H, which is the physical energy operator. To avoid confusion we prefer to 
speak of low or high momentum states. 

As a unitary transformation leaves eigenvalues invariant it is quite natural to perform a suitable 
unitary transformation of the uncorrected many-body states |0[v]) to another basis |0[v]), which 
already includes correlations to a certain extent so that the Hamiltonian matrix becomes more 
band-diagonal and is not scattering to high momentum states anymore. Whenever many-body 
states include short-range correlations we mark them by ~ like the fully correlated eigenstates 
of the Hamiltonian H which includes short-range repulsive and tensor potentials. 

In this contribution we discuss in depth the Unitary Correlation Operator Method (UCOM) 
which tries to achieve this goal by an explicitly given unitary operator that transforms the Hamil- 
tonian and all other observables to effective operators that include the effects of the short-range 
correlations. 

Before discussing effective interactions in general and UCOM in particular, we would like to 
address a few more general issues related to the nuclear potential and effective interactions. 

1.1. Many-body potentials 

In recent years it has become clear that the nuclear many-body problem needs at least three- 
body forces to reproduce data with sufficient precision 

A A A 

H = T"UvUi + vUi^ + --- = ;^T,+ j;V,,+ Vo. + .... (5) 

;=1 y>'=l k>j>i=\ 

An irreducible three-body potential V|^^j^ is the remaining part of the Hamiltonian acting in three- 
body space that cannot be described by a sum over one-body kinetic energies and two-body poten- 
tials. Likewise, irreducible four-body potentials are the remaining parts of the Hamiltonian acting 
in four-body space that cannot be described by a sum over two-body and three-body potentials and 
so on. 

A typical example of an irreducible or genuine three-body potential is the Fujita-Miyazawa 
interaction [jvl], which arises from the restriction to proton and neutron as elementary degrees of 
freedom. Due to the strong coupling with the pion, a nucleon can be converted into a A-resonance, 
which then has to decay back into a nucleon, because the energy of the total system is not sufficient 
to produce a A on its mass shell. This intermediate virtual excitation is not contained in the 
nucleonic Hilbert space, but is effectively already included in a two-body potential that is fitted to 
data. The left-hand graph displayed in Fig.[Tl where the third nucleon is only a spectator, therefore, 
is a contribution that makes the two-body potential an effective one with respect to excluded A- 
degrees of freedom from the nucleonic Hilbert space. In the right-hand graph the pion is coupling 
to a third nucleon. Thus this interaction cannot be written as an effective two-body potential. It 
has to be seen as a genuine three-body potential. 

Already from Fig. [T] it is obvious that two-body and three-body forces cannot be treated inde- 
pendently as they may originate from the same physical process. Furthermore, this simple example 
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Figure 1 : A-resonance as intermediate off-shell excitation in a three-body system. 



shows that omitting parts of the many-body Hilbert space that are in principle reachable with the 
strong interaction leads to effective operators and induces many -body forces. This effect will also 
occur, and is of major concern, in the following sections when the nucleonic Hilbert space itself is 
truncated to low-energy states. 

1.2. Off-shell properties 

The elastic nucleon-nucleon scattering cross section is represented in terms of phase shifts 
and mixing angles as a function of energy. As cross sections are measured far away from the 
interaction region, these quantities give only indirect information about the nuclear forces and 
cannot determine the potential in a unique way. In a stationary scattering state the two nucleons 
have a sharp energy and the modulus of the asymptotic relative momentum q of the ingoing 
wave has to come back on-shell to the same momentum q in the outgoing wave, both related to 
the eigenenergy by Eq = q^/(2/j). 

In an interacting many-body system the total energy E„ corresponding to an eigenstate of the 
Hamiltonian is sharp but a pair of particles within the system has neither a sharp relative energy 
nor a sharp momentum. While they feel their mutual interaction they are also interacting with 
other particles exchanging energy and momentum. Therefore, the interaction is tested for different 
in- and outgoing momenta q and q' that are not connected by any on-shell relation. Studying 
many-body systems will, therefore, give information on the off-shell behavior of the nuclear force. 

As we will see in Sec. [2] a unitary transformation that acts only at short distances and does 
not affect the asymptotic behavior, like the UCOM correlator, leaves the phase shifts untouched 
but creates a new potential that is not distinguishable from the original one by measurements of 
the elastic cross section. As one can devise an infinite number of such transformations, there are 
infinitely many phase-shift equivalent potentials. These, however, differ in their off-shell behavior. 
Therefore, scattering data will never be able to determine the nuclear force uniquely. In Sec.|4]and 
[5]we show that phase-shift equivalent potentials with different off-shell properties can produce very 
different results when applied in many-body systems. One of them produces with increasing mass 
number a dramatic overbinding in a variational calculation, where the calculated binding energy 
is, anyhow, only an upper limit. In order to be sure that the limited Hilbert space used is not 
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the reason that the other potentials are not overbinding as well one needs many-body approaches 
models that can solve the many-body problem exactly for a given Hamiltonian. 

One goal of modern nuclear structure theory is to disentangle as much as possible the off- 
shell nature of the potential from effects originating in the use of restricted many-body Hilbert 
spaces. However, one should always keep in mind that in any theory, not only in nuclear structure, 
any interaction is an effective low-momentum interaction in the sense that it is constructed for 
the degrees of freedom and Hilbert space one is using. With increasing momentum transfer one 
always opens up new degrees of freedom not contained in the original model space. 

Exact benchmark calculations for the three- and four-body system dsLls, 1^ 11, 12] were very 
helpful in this respect. They also showed that three-body forces are needed, though different 
ones for different phase- shift equivalent two-body potentials. For example the contribution to the 
binding energy of "^He coming from three-body terms complementing the Argonne VI 8 potential 



are about 50% larger than those accompanying the CD Bonn potential al3\\ . see Fig. |2T1 The CD 
Bonn potential has a softer short-range repulsion than the AVI 8 but contains a radial momentum 
dependence which is absent in Argonne VI 8. Glockle and collaborators have shown in a quite 
general way that off-shell properties of two-body interactions can be traded against three-body 
forces {T^ . We will encounter this effect again when we perform the UCOM and SRG similarity 
transformations of the Hamiltonian. 

1.3. Effective potentials 

Besides the conceptual problems of deriving and defining the nuclear interaction there is also 
the already mentioned technical problem that the solution of the many-body Schrodinger equa- 
tion (HI) requires a representation of the many-body eigenstates. Since realistic nuclear interac- 
tions induce various kinds of correlations, in particular short-range correlations, one possibility is 
to work with many-body states that can represent the correlations. In that case the main numerical 
effort goes into calculating the matrix elements of the Hamiltonian. Exarnples are the exact Fad- 



deev and Faddeev-Yakubovsky equations for the 3- and 4-body system lllSn . the hyperspherical 



ivpersp 

harmonics basis 1160, or the Green's Function Monte Carlo method (GFMC) [ll7Lll8L]l9[ EoL l21ll. 
The other possibility is to represent the eigenstates of the Hamiltonian with many-body basis 
states |0[y]) which are chosen such that the numerical effort for calculating the matrix elements 
<0[v]| T"] |0[^]), <0[,]| V^J^ |0^]), (O^l Vj^^jljj^ |0^]) of one-, two-, and three-body operators is min- 
imized. In that respect the best choice are antisymmetrized product states (Slater determinants). 
However, the strong short-range correlations induced by realistic nuclear forces cannot be repre- 
sented in a reasonable way by a product state basis. For example the four-body system "^He, when 
represented in harmonic oscillator states, needs a three-body Hilbert space of dimension more than 
10^ to get a converged ground-state energy with the Argonne VI 8 potential, see Sec. HI Therefore, 
so called effective interactions for truncated Hilbert spaces are introduced. At present there is no 
other possibility for mass numbers larger than about 12. 

1.4. Unitary approaches 

In a product basis representation the nucleon-nucleon interaction scatters to energetically very 
high-lying basis states. Or in other words, the short-range repulsive and tensor correlations imply 
components in the many-body state with large relative momenta, which necessitate very large 
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many-body Hilbert spaces in order to accommodate the correlations in this basis. To still work with 
the numerically convenient product basis one uses so called effective interactions. The effective 
interaction should decouple the Hilbert space containing high momenta from the one with low 
momenta, which can be represented more easily by product states. 

The most straightforward and intuitive way is to consider a unitary transformation of the prod- 
uct basis |0[y]) _ 

|Om) = U |Om) (6) 

to render the Hamiltonian matrix into a more diagonal form not connecting low and high momenta. 
The new basis \^[v]) should contain already the main properties of the short-range correlations like 
a depletion of the many-body wave function 

<^l, ^2, • • • , ^a|Om> with = (4, (Tt, T,-) (7) 

whenever the distance \xi - Xj\ between two nucleons of the many-body system is within the range 
of the repulsive core of the interaction. Similarly their spins (T/, (Tj should be aligned with fij = 
Xj - Xj if their isospin is in the T = configuration, because this correlation produces binding 
energy, see Fig. [21 

This unitary transformation can be used to define an effective Hamiltonian through the simi- 
larity transformation 

Heff =U"^HU = U^HU. (8) 

The unitary transformation fulfills the principal requirement formulated in Eq. ([3]) that the eigen- 
values for the energy are invariant. Later approximations should be checked regarding this aspect. 
Another advantage is that orthonormality relations between eigenstates of H or between basis 
states are not changed 

^ (^J^,«) = (^«|U-^U|T,„) = <¥J»PJ = 5„,,„ 
(OmIO^,]) = (OmI U-^U |0[^]) = (OmIO^,]) = . (9) 

Any unitary transformation can be written in terms of a hermitian generator G as 

U = e-'^ . (10) 

If G = T/j=i Sj is a one-body operator no correlations are induced because in that case U = 
YVj=\ just transforms each single-particle state in |<l>[v]) independently. Therefore, the gen- 
erator G has to be at least a two- body op erator. This implies that the unitary transformation of a 



one-body operator T^'^ yields [1221 l23l l24h 



Ut u = + Tf^i + Tf^i + . . . . (11) 



The transformation of a two-body operator V'^^ 



NN 



yields a new two-body operator V[^^ and additional operators of higher-order. Likewise a trans- 
formed three-body operator results in a new irreducible three-body operator with its higher-body 
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companions. The above the operators indicates that effects from short-range correlations are 
now moved to the operators in the sense that the potential Vf^^ is softer or less repulsive at short 
distances than the original V|^^. This "taming" of the potential comes at the expense of introducing 
two-, three- and more-body terms originating from the kinetic energy as well as from the potential. 
The different contributions will be discussed and shown in Sec. [21 

An important message that holds for any method of deriving effective interactions is that the 
effective Hamiltonian contains irreducible n-body interactions, where n goes in principle from 1 
to A. For example, if we start from a Hamiltonian with two- and three-body forces like in Eq. ([5]), 
we obtain 

U-^ HU = Tl" + (t^ + N'^^) + (Tl^l + V|,^] + V|,^] ^) + . . . . (13) 

The hope is to keep the n> 2 body terms small in order to reduce the numerical effort. 

One important advantage of formulating the effective Hamiltonian through a unitary similarity 
transformation is that any other observable B can and should be transformed the same way. An 
arbitrary matrix element of B between two eigenstates of the Hamiltonian can be written as 

<^„| B |¥,„) = <¥„| U'^B U |¥,„) = <^„| Beff |¥,„) (14) 



which implies the definition of Bgff as 

Beff = U^BU. (15) 

Again Bgff becomes a many -body operator even if B is an one-body operator, see Eq. (fTTI) . In 
section[2]and|3]approximations in the calculations of Hgff are introduced. The same approximations 
can also be applied to calculate Bgff . 

1.5. Jastrow ansatz 

An early attempt to incorporate short-range correlations in the many-body state was proposed 



by Jastrow [|25L |26|] 



|0, = <5 f] &U C?;, Tu Tj) |0) (16) 

where in modem applications the Jastrow correlation functions / contain besides the relative dis- 
tances ?ij = - ij, also operators depending on spins and isospins to account for correlations 
other than the short-range repulsion, like the tensorial ones. For distances \rij\ much larger than 
the range of the interaction /(/^;, . . . ) approaches one. As the symmetrized product runs over all 
particle pairs, the transformation from the uncorrelated state |0) to the correlated |0, J) is an A- 
body operator, like in all other approaches. Even three body products /(x,, Xy, x^, . . .) are being 



used[|i7L|20L|26|] 



Although very intuitive, the Jastrow ansatz implies a huge numerical effort in calculating ex- 
pectation values and matrix elements. Even the norm of |0, /) cannot be calculated analytically. 
As the transformation is not unitary, configuration mixing calculations encounter even greater nu- 
merical challenges. 
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1.6. Projective approaches 

In the past and still nowadays many effective interactions are based on projection methods. In 
this case one is looking for an effective Hamiltonian that should have the same eigenvalues as 
the full Hamiltonian H and its eigenstates should be the projections of the exact eigenstates I*?,,) 
onto the model space: 

PHeff |0„) = E„ |(D„) and |0„) = P |¥„) . (17) 

To simplify the equations one writes the Hamiltonian H = Hq + V and the effective Hamiltonian 
Hgff = Ho + Vgff as a sum of Hq and a residual interaction V and Vgff, respectively. The formal 
solution for Hgff that fulfills the requirement (fTTl) is actually energy-dependent and given by 

Vgff(£) = V + V ^ Vgff(£) , (18) 
{E - Ho) 

where the energy E has to be taken as the unknown eigenvalue E„ of the full problem. This means 
that each eigenstate in principle corresponds to a different effective Hamiltonian. As the difference 
between E and eigenvalues of Hq that belong to the Q-space enters Eq. (fTSl) . the energy dependence 
will be weak if the lowest Q-space eigenvalue is far away from the considered £■„. Note that Veff(£') 
in this equation is not restricted to the P-space but connects P- and Q-space. However, the effective 
interaction is actually PVgft(£')P. The full Vgff(£') is needed when one calculates other effective 
observables (see below). Again Eq. (fT8l) yields an effective Hamiltonian that contains many -body 
operators. 

When solving this equation in two-body space one obtains the Brueckner G-matrix 

G{oj) = V + V— ^ G{aS) , (19) 
(oj - Ho) 

where the choice of the starting energy a» and the Pauli projection operator q acting in two-body 
space is not uniquely defined. After having solved the G-matrix equation (fT^l ) one uses G{co) in 
the A-body space, which constitutes a low density approximation. 

Although Eq. (fTSl) appears to be suited for a perturbation expansion, this is not possible if the 
strong short-range correlations of the nuclear interaction have to be renormalized. At least infinite 
partial summations are needed. The Brueckner G-matrix is such a partial summation, namely over 
the so called ladder diagrams. 

Another task (often not considered) is that not only the Hamiltonian should be transformed to 
an effective one, but also all other observables B such that 



(0„|Bg,|0.) <^JPBg,-P|^J ^^^^^1^1^^^ ^^^^ 



<(D„|0„)i/2<(D„,|0„,)i/2 <0„|0„)i/2<0,JcD„,) 



1/2 



for the low lying many -body states of interest. This could in principle be done by using the A-body 
Bethe-Goldstone equation 

l^«> = ( 1 + 7F-^TFT Veff(£„)) , (21) 

^ (En - Ho) ' 
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which also follows from Eq. (fTT]) . It reconstructs the correlated many-body state |*P„) from its 
projected part |0„) = P When inserting Eq. (|2T]) into Eq. (I20l) to obtain the effective op- 
erator one does not need much imagination to see that this is a rather impracticable method and 
susceptible to errors in approximations of the matrix elements of YefiiEn) that connect P- and Q- 
space. Within the G-matrix approximation one uses the analoguous Bethe-Goldstone equation, 
where YesiE) is replaced by G(oj) and Q by q. The G-matrix method is not discussed further in 
this contribution. 

Unitary approaches are more transparent and not energy dependent. A unitary approach is 
more robust as there is some freedom in the choice of U, because only a few low lying eigenvalues 
and eigenstates of PHgffP = PU 'H UP are needed. A not so optimal choice of U can be compen- 
sated by a larger P-space. Unitary approaches can even generate effective interactions that are by 
construction phase-shift equivalent with the original interaction. 

1.7. Plan of the review 

In Section |2l the concept of the Unitary Correlation Operator Method (UCOM) is laid out. 
Central and tensor correlators are introduced and their properties are illustrated with descriptive 
examples. Section [3] discusses the relation to the Similarity Renormalization Group (SRG) ap- 
proach, which is more general, but in its application to nuclear effective interactions it is based on 
the same physical assumptions. From SRG we deduce a set of optimal correlation functions for 
UCOM and compare the resulting effective low-momentum interactions. Section |4] steps into the 
many-body Hilbert space using the No-Core Shell Model (NCSM) to probe the effective interac- 
tions developed in the preceding sections. Properties, like masses, spectra, radii, magnetic dipole 
and quadrupole moments, are investigated for nuclei with mass numbers A = 3 -7. In Section[5]we 
move on to nuclei along the whole nuclear chart by means of the Hartree-Fock approximation and 
low-order many-body perturbation theory. We show that the UCOM concept is very versatile not 
only for the understanding of the nuclear many -body system but provides a general method to treat 
short-range correlations in a quantitative way. We also discuss evidence for missing three-body 
forces or density dependences, which seem to be needed to obtain the correct nuclear saturation 
properties for large systems. Finally we close with summarizing conclusions. 
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Figure 2: Argonne V18 potential in the S - l,T - channel. The lower curve is obtained when the spins are 
aligned parallel to the distance vector of two nucleons whereas the upper curve is obtained when the spins are oriented 
perpendicular to the distance vector. This spin dependence of the potential is caused by the tensor force. The relative 
momentum q is chosen to be zero, therefore, the spin-orbit force does not contribute. 



2. Unitary Correlation Operator Method (UCOM) 

As discussed in the introduction, short-range central and tensor correlations pose a major chal- 
lenge for the solution of the nuclear many-body problem. These correlations are induced by cor- 
responding features of realistic nucleon-nucleon forces that are illustrated in Fig. [21 where the 
Argonne VI 8 potential is plotted as a function of the distance of the nucleons and the orientation 
of the spins relative to the distance vector. At short distances, r < 0.5 fm, the potential is strongly 
repulsive, it has a repulsive "core" which induces short-range central correlations. At larger dis- 
tances the potential shows a pronounced dependence on the orientation of the spins. If the spins 
are aligned perpendicular to the distance vector the potential is almost flat, whereas the potential 
is attractive with a minimum at r « 1 .0 fm if the spins are aligned parallel to the distance vector. 
This difference is caused by the tensor force, which originates mainly from the one-pion exchange 
part of the potential. 



In the Unitary Correlation Operator Method 11271 l24i |23| . 12211 these correlations are imprinted 



into a many-body state by means of a unitary correlation operator C: 

1^) = C 1^) . (22) 

Here the uncorrelated many-body state is a "simple" state that cannot represent the short-range 
correlations. This can be a Hartree-Fock Slater determinant or a basis state of the No-Core Shell 
Model. 

Short-range correlations are important when the particles come close. If the density is low 
enough, the probability for finding three nucleons within the correlation volume defined by the 
correlation range will be small. In that case the assumption that the short-range correlations are of 
two-body nature will be a good approximation. At nuclear saturation density po = O.lVfm"^ the 
mean distance between nucleons is 1.8 fm while the repulsive core sets in around 0.5 fm. There- 
fore, we assume that the correlations are essentially independent of the environment and we use 
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a state-independent ansatz for the correlation operator C. This implies that the short-range cor- 
relations for low-lying states are very similar in all nuclei. The unitary correlation operator C 
describing this transformation is given in an explicit operator form, independent of a particular 
representation or model space. The correlation operator is optimized for the lowest orbital angular 
momentum L in each spin-isospin channel but it does not explicitly depend on L. The UCOM 
approach can, therefore, be also used in many-body approaches that do not use basis states of 
good angular momentum. Alternatively, one could define different correlation operators for each 



partial- wave channel. In the S -wave UCOM approach proposed by Myo and Toki 112 8h only central 
correlations in L = channels are considered. 

When calculating expectation values or matrix elements of some operator A we can either 
evaluate the bare operator A in the correlated states I^F) or we can use a correlated operator A 
defined through a similarity transformation 

A = C"'AC = CUC (23) 

with the uncorrected states. 

<¥| A = <¥| C'^AC = <*F| A . (24) 

Due to the unitarity of C, the notions of correlated states and correlated operators are equivalent 
and we may choose the form that is technically more advantageous. 

In the case of the nuclear many -body problem, the unitary correlation operator C has to account 
for short-range central and tensor correlations as explained above. In the UCOM approach we 
explicitly disentangle these diiferent types of correlations and define the correlation operator as a 
product of two unitary operators, 

C = CnC, , (25) 

where Cn describes short-range tensor correlations and C,- central correlations. Each of these 
unitary operators is expressed with hermitian two-body generators 

Cn = exp[ - i ^ gn,,;] , C, = exp[ - i ^ g,,,y] . (26) 

The details of the generators g,. and gn will depend on the particular nucleon-nucleon interaction 
under consideration. 



2.1. Central correlations 

The central correlations are induced by the repulsive core of the central part of the interaction 
which tries to keep the nucleons apart from each other. The two-body density in the correlated 
many-body state will be strongly suppressed at short interparticle distances, i.e. in the range of 
the repulsive core, and it will be enhanced at larger distances, where the potential is attractive. 
This can be achieved by a distance-dependent shift in the relative wave function for each pair of 
nucleons. The generator g,. is constructed such that it performs these shifts in a unitary way. The 
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shifts are generated by the projection of the relative momentum q = ^[pi - P2] onto the distance 
vector = xi - X2 of two nucleons: 



with r = The amplitude of the shift — large shifts at small distances within the core, small or no 
shifts outside the core — is described by a function ssrir) for each spin-isospin channel. The de- 
tailed form of the function ssrir) will depend on the potential. Its determination will be discussed 
in detail in Sec. 12.51 The full generator for the central correlations is written in a hermitized form 
as ^ 

gr = ^2^^r ssAr) + ssrir) qA Ust , (28) 

S,T 

where Ust is the projection o per ator onto two-body spin S and isospin T. Similar generators were 



already used by Ristig et al. ||29L |30L |3l|] for the description of the central correlations induced by 
hard-core potentials. 

2.2. Tensor correlations 

The correlations induced by the tensor force entangle the orientation of the spins with the spa- 
tial distribution of the nucleons to optimize the contribution of the tensor force. We construct the 
tensor correlation operator in such a way that it will only act on the orbital part of the relative wave 
function of two nucleons. This can be achieved by using only the "orbital part" qn of the relative 
momentum operator, which is obtained by subtracting the radial part of the relative momentum 
operator (used in the generator for the central correlations) from the full relative momentum oper- 
ator 

qn = q - - • q,- = ;^(L X ?- ?x L) (29) 
r 2r'^ 

with L = i^x ^. Like the tensor operator Si2(;, ;) in the tensor force, the generator for the tensor 
correlations is the scalar product of an operator of rank 2 in coordinate space (constructed from the 
relative distance vector r and the "orbital momentum" qn) with the operator of rank 2 in spin-space 



i23n . As for the central correlator the amplitude of the shifts will depend on the distance of the 
nucleons and the potential under consideration. The amplitude is given by the tensor correlation 
function ??7-(r). In cartesian notation we can write the full generator as 

ga = J]^T(r)SnirAn)U,T (30) 

T 

using the general definition for a tensor operator of rank 2 

Si2(a, b) = |[((?i • ^(c?2- b) + ((?i ■ b)((?2- ^] - K^i •^2)(a ■ b + b ■ ^ . (31) 
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2.3. Correlated wave functions 

As explained above, because the correlation operator acts only on the relative motion of a nu- 
cleon pair, the center-of-mass motion of the nucleon pair is not affected. When discussing corre- 
lated two-body wave functions we can, therefore, restrict ourselves to the relative wave functions. 
For the uncorrected relative wave function we assume LS -coupled angular momentum eigenstates 
\(f)(LS)JM TMj) with the radial wave function (p{r). The correlation operators do not depend on 
M and Mj and we will omit these quantum numbers in the following. 

The central correlator = exp(-i g^) affects only the radial part of the wave function and 
leaves the orbital part of the wave function and the spin and isospin unchanged. In coordinate rep- 
resentation the correlated wave function can be rewritten as a norm-conserving coordinate trans- 
formation [|221 



{r{L'S)JT\ C- \(P{LS)JT) = ^JR^ (l>{R-{r))5 u i 



r 



(32) 



{r{L'S)JT\Cl \(P{LS)JT) = ^JWJr) ((>{R^{r))6L' 



where R+{r) and R-{r) are mutually inverse, R+(R+{r)) = r. The correlation functions R+(r) and 
R-{r) are related to the function s(r) used in the generator (|28l) through the integral equation 

= ±1. (33) 



For illustrative purposes the correlation functions can be approximated as R±(r) « r + s{r). 

In the LS -coupled basis the application of the tensor correlator Cn can be expressed easily. 
The tensor operator 812(1% qn) used in the generator has only off-diagonal matrix elements in the 
LS -coupled basis 

<(7 ± 1, l)7r| Si2(?, qn) 1(7+1, 1)7^) = ±3i /7(7TT) . (34) 

Within a subspace of fixed 7 one can, therefore, calculate the matrix exponential and thus the 
matrix elements of the full tensor correlator Cn- 

The tensor correlation operator will have no effect for states with L = J, whereas states with 
L = 7 + 1 will be connected to states with L = 7 + 1 . The strength of the mixing is governed by 
&(r) 



{r(L'S)JT\Cn \((>iLS)JT) = 



cos 6j(r)(p(r) 
± sin ej{r)(p{r) 



U =L = J 

U =L = J +\ 

L' = 7 + 1,L = 7 + 1 



(35) 



where we use the abbreviation 



ej{r) = 3 ^[J{JTT) &(r) . (36) 

Combining central and tensor correlations we end up with the following expression for the 
fully correlated wave function in coordinate space: 

'^,JWJr)cf>(R_(r)) ■L'=L = J 



{r(L'S)JT\ CnC, \(P{LS)JT) = \ cos OAr)^ ^JWJr) (f>(R-(r)) 

■.smejir)^^/WJrj(f>iR.ir)) 



U =L = J +1 

U = J +l,L = J+\ 



(37) 
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Figure 3: Central and tensor correlations in the deuteron channel. Starting from an uncorrelated trial wave function (a), 
first the central (b) and then also the tensor correlated wave functions (c) are shown together with the corresponding 
central (d) and tensor correlation functions (e). 



The wave function in momentum space can be obtained by Fourier transformation of the coor- 
dinate space wave function (|37]) 

{q(L'S )JT\ CnC, |0(L5 )JT) = ^i^' J drr^ iu{.qr){r{L' S )JT\ CnC, |0(L5 )JT) , (38) 

where we use momentum eigenstates normalized as {q{LS)JT\q'{L'S)JT) = ^5{q - q')6i>i. 

The deuteron wave function provides an illuminating example for the role of central and tensor 
correlations. We will start from an uncorrelated state \^o(LS)JT) = |(/io(01)10), which is a pure 
S -wave state with the spin-isospin quantum numbers of the deuteron. The radial wave function 
(po(r) shall not contain short-range correlations induced by the repulsive core. Figure |3] shows the 
uncorrelated L = radial wave function. Applying the central correlator C,- with the correlation 
function R+(r) leads to a wave function with a correlation hole at short interparticle distances. The 
application of the tensor correlation operator Cn in a second step generates the Z)-wave compo- 
nent in the wave function, which depends on the tensor correlation function {^(r). If we assume 
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an uncorrelated wave function that is purely S -wave, the entire D-wave component has to be gen- 
erated by the tensor correlator. This can be achieved using a very long-ranged tensor correlation 
function (as given by the dashed curves in Fig. [3]). But this is not the idea of the unitary correlation 
operator method, where we only want to describe the state-independent short-range correlations 
by means of the correlation operators. Long-range correlations, like those governing the outer 
part of the deuteron wave function, should be described explicitly by the many-body approach. 
The solid curves show the Z)-wave component obtained with such a short-range tensor correlator 
(7^ = 0.09 fm^ as discussed in Sec. 12.51) . 



2.4. Correlated operators 

Correlated Hamiltonian - Central Correlations 

In the two-body system the unitary transformation with the central correlator C,- can be ex- 
pressed analytically for any operator that can be written as a function of relative distance i^ and 
relative momentum ^. 

The most important example is a Hamiltonian consisting of kinetic energy and a realistic 
nucleon-nucleon interaction given in a generic operator form 

^ = Z ^[^i'WO,' + OpVpij)] (39) 
p 

with 

Op = {1, ((?i-(?2), q', q'((?i-(?2), L', LVrf?2), ^^^^ 
(L-^), Si2(J,J), Si2(L,L)}®{l, (frr2)}. 

Here we only consider non-local terms up to quadratic momentum dependence. Such terms appear, 



e.g., in the Nijmegen 113211 or operator representations of the Bonn A/B potentials A quadratic 
momentum dependence of the form q^ can be expressed by the q^ and 1? terms contained in (l39l) . 
Charge dependent terms in the interaction are not explicitly discussed, but are included in the 
calculation of matrix elements and the many-body calculations based on them. 

For the formulation of the correlated Hamiltonian in two-body space, we start with an initial 
Hamiltonian given by 



H = Hint + Ten, = Tint + V + Tc,„ = T, + Tfj + V + Ten, , 



(41) 



where we have decomposed the kinetic energy operator T into a center of mass contribution Ten 
and an intrinsic contribution Tint which in turn is written as a sum of a radial and an angular part 



1 , _ 1 p 



Applying the central correlator C,- in two-body space leads to a correlated Hamiltonian consisting 
of the one-body kinetic energy T and two-body contributions for the correlated radial and angular 
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relative kinetic energy, TJ. and , respectively, as well as the correlated two-body interaction 

Cj: H C, = T + T|^] + Tg^ + yf^^ . (43) 

The explicit operator form of the correlated terms can be derived from a few basic identities. The 
unitary transformation for the relative distance operator r results in the operator-valued function 

ClrCr = R^{r). (44) 
Because of the unitarity of the correlation operator C,. an arbitrary function of r transforms as 

/(r) Cr = f(Cl r C) = f(R,(r)) . (45) 

The interpretation of the unitary transformation in terms of a norm-conserving coordinate trans- 
formation r R+(f) is evident. For the radial momentum operator one finds the following 
correlated form ll22ll 

Cl q, C, = —L= q,. -jL= . (46) 
With this we can express which enters the radial part of the relative kinetic energy as 

C;q?C = i(,J^.^,J) + TO, (47) 

which consists of a transformed momentum dependence plus an additional local term depending 
only on the correlation function R+{r) 

lR'!{rf R'!'{r) 

W(r) = — — i-— . (48) 

^ ' 4i?;(r)4 2i?;(r)3 ^ ' 

All other basic operators, such as L^, (L • S), Si2(;, \) commute with the correlation operator C,- 
and are, therefore, unchanged by the central correlations. 

Based on these elementary relations we can explicitly construct the two-body contributions to 
the correlated kinetic energy. For the radial part we obtain using (1471) 

T- = C;T,a - T, = . ^qj) . i-H'(r) (49) 

with a distance-dependent effective mass term 

1 _ 1 / 1 X 

The two-body contribution to the correlated angular part of the kinetic energy involves only the 
basic relation (|46|) and gives 

Tg] = CjTnC,- - T^, = ,r\-^ (51) 

2yun(r) r2 
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with a distance-dependent angular effective mass term 

Ifiair) 2iu^R+(ry ' 

The momentum dependent terms of the nucleon-nucleon interaction (l39l) transform in a similar 
manner like the kinetic energy. Using (l46l) and (l47l) we obtain 

C; i(q;V(r)+V(r)qJ) C = TO-V'(^,(r))||l . (53) 

For all other terms of the NN-interaction (l39l) the operators Op commute with the generator and 
we only have to transform the radial dependencies 

Cj y(r)OpC, = y(i?+(r))Op. (54) 

Many of the other relevant operators, e.g. the quadratic radius or transition operators, can be 
transformed just as easily. 

Correlated Hamiltonian - Tensor Correlations 

The transformation of the Hamiltonian with the tensor correlation operator Cq is more in- 
volved. In general, it can be evaluated via the Baker-Campbell-Hausdorff expansion 

CiACn = A + i[gn,A] + ^[gn,[gn,A]] + --- . (55) 

In some cases the series expansion will terminate after a finite number of terms. A trivial case is 
the distance operator r, which commutes with the tensor generator gn and is thus unchanged by 
the transformation 

C|,rCn=r. (56) 

For the radial momentum operator q^, the expansion (l55l) terminates after the first order commuta- 
tors and we obtain the simple expression 

qr Cn = q. - ^'(r) Si2(?, qn) ■ (57) 

For the tensor correlated quadratic radial momentum operator the series terminates after the first 
two terms and we obtain 

Cn = q? - [&'(r) q, + q, S^i?, 4n) + [^'(r) 812(1^, , (58) 

where S\2{?,qnf = 9[§^ + 3(L • ;?) -1- (L • By applying the tensor correlator to the kinetic 
energy we have generated momentum-dependent tensor operators as well as "conventional" spin- 
orbit and tensor terms (from the (L • term). The correlated kinetic energy is no longer a central 
operator. 
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For all other operators of the interaction (l39l) that depend on angular momentum, the Baker- 
Campbell-Hausdorff series does not terminate. Through the commutators additional tensor opera- 
tors are generated. At first order we obtain 

[gii, Snil f)] = i'?(r)[-24ni - 18 (L ■ ^) + 3 Suil ?)] 
[gn, (L • S)] = i7?(r)[-Si2(qa, qn)] ^^^^ 

[gi2,Si2(L,L)] =ii?(r)[7Si2(qn,qn)] , 
where we use the abbreviation 

Sn(4n, 4n) = 2r%2(4n, 4n) + SniU L) - 5 Si2(?, ?) ■ (60) 

The next order generates terms of higher order in orbital angular momentum, e.g., a L^(L • §) 
term. In practice we have to truncate the Baker-Campbell-Hausdorff expansion to some finite set 



of operators [|24l1 . In principle the contributions of the higher-order operators will become more 
important with increasing angular momenta. On the other hand the contributions of the correlated 
Hamiltonian, which are of short range, will be overwhelmed by the centrifugal barrier from the 
one-body kinetic energy for large angular momenta. 

Note that matrix elements of the correlated Hamiltonian using angular momentum eigenstates 
are calculated by applying the tensor correlation operator onto the basis states (see Sec. 12.31) . which 
does not require approximations. 

2.5. Optimal correlation functions 

The central and tensor correlators depend on the correlation functions s(r) and ??(r) in the dif- 
ferent spin- and isospin-channels. We now have to determine these correlation functions for a 
given nuleon-nucleon potential. One important question, that was already raised in the discussion 
of the deuteron wave function, is the separation between state-independent short-range correla- 
tions, which we want to describe by the correlation operator, and long-range correlations that 
should be described explicitly by the many-body approach. 

The most convenient procedure to determine the correlation functions is based on an energy 



minimization in the two-body system 1271 1231] . For each combination of spin S and isospin T we 
compute the expectation value of the correlated energy using a trial state with the lowest possible 
orbital angular momentum L. The uncorrelated radial wave function should not contain any of 
the short-range correlations, i.e., it should resemble the short-range behavior of a non-interacting 
system. In the following we will use a free zero-energy scattering solution (piir) oc r^. Other un- 
correlated trial wave functions, e.g. harmonic oscillator eigenfunctions, give very similar results. 

For practical reasons the correlation functions are represented by parametrizations with typ- 
ically three variational parameters. The drop-ofi" can be well-described by a double-exponential 
decay with variable range. For the short-range behavior, several difi'erent parametrizations have 
been compared. For the Argonne VI 8 potential, the following two parametrizations for the central 
correlation functions have proven appropriate: 

Rl(r) = r + a (r/fif exp[- exp(r/yS)] , 

IT (61) 
R'^(r) = r + a[l - exp(-r/7)] exp[- exp(r/yS)] . 
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T 


Param. 


a [fm] 


ySLfm] 


7 [fm] 


n 








II 


0.7971 


1.2638 


0.4621 







1 


I 


1.3793 


0.8853 




0.3724 


1 





I 


1.3265 


0.8342 




0.4471 


1 


1 


II 


0.5665 


1.3888 


0.1786 





Table 1 : Parameters of the central correlation functions R+{r) specified in ( |6TI ) for the Argonne V 1 8 potential obtained 
from two-body energy minimization. 



Which of these parametrizations is best suited for a particular channel will be decided on the basis 
of the minimal energy alone. All parametrizations allow only outward shifts by construction. This 
is different from the correlation functions that are obtained by the SRG mapping procedure (see 
Sec. l3.6|) . For the tensor correlation functions the following parametrization is used 

i^(r) = a [1 - exp(-r/7)] exp[- exp(r/;0)] . (62) 

The 5=0 channels are only affected by the central correlators. Their parameters are de- 
termined by minimizing the energy for the lowest possible orbital angular momentum state, i.e., 
L = 1 for r = and L = for T = 1 , respectively, 

£oo = <0i(lO)lO|CjH,„tC,|0i(lO)lO) , 
£oi = <0o(OO)Ol|C:H,„tC,|<^o(OO)Ol) . 

For S = 0, r = 1 the minimization of Eq\ by variation of the parameters of the central correlation 
function is straightforward. The resulting parameters are summarized in Tabled] For S = 0, T = 
the potential is purely repulsive and, therefore, the energy minimization leads to central correlation 
functions of very long range. In order to avoid this pathology we employ a constraint on the 
strength of the correlation function given by 

= J drrHRAr)-r). (64) 

The value of this constraint on the central correlation function for the S = 0, T = channel is 
fixed to Ir^ = 0.1 fm'* in accordance with typical values in other channels. 

For S = 1 we also have to consider tensor correlations and the parameters of the central and the 
tensor correlation functions have to determined simultaneously. For T = the energy is defined 
by the matrix element with L = states 

£io = <<^o(01)10|C+C;,Hi„tCnC,|<^o(01)10) . (65) 

In the r = 1 channel the lowest possible orbital angular momentum is L = 1. From angular 
momentum coupling we obtain 0, 1, and 2 as possible values for J. Therefore, we define the 
energy functional used for the minimization procedure by averaging over the angles, which is the 
sum over all three possibilities with relative weights given by 27 -I- 1 

En = i<0i(l 1)011 Cj:Hi„tC,|<^i(ll)01) 

+ |<0i(l 1)1 l|CjHi„tC, 101(11)11) (66) 
+ |<0i(l 1)211 CjC^ Hint CnC, 101(1 1)21) . 
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Figure 4; Optimal central correlation functions R+{r) — r for the Argonne V18 potential according to the parameters 

given in Tab. [T] The curves correspond to the different spin-isospin channels; S - 0,T - 1( -), 5 - 1,T = 

( ), 5 = 0, r = ( ), and 5 = 1, r = 1 ( ). 
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Table 2; Parameters of the tensor correlation functions defined in ( l62T i for the Argonne V 1 8 potential with different 
values I§ for the range constraint obtained from two-body energy minimization. 
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Figure 5: Optimal tensor correlation functions §{r) for different values of the range constraint (a) Correlation 

functions for T = with h = 0.06 ( ), 0.09 ( ), and 0.15 fm^ ( ). (b) Correlation functions for T = 1 

with h = -0.01 ( ), -0.03 ( ), and -0.09 fm^ ( ). 



As mentioned earlier, the long-range character of the tensor force leads to long-range tensor 
correlations. However, long-range tensor correlation functions are not desirable for several rea- 
sons: (0 The optimal long-range behavior would be strongly state-dependent. Therefore, our goal 
of extracting the state-independent, universal correlations forbids long-range correlation functions. 
(//) The two-body approximation would not be applicable for long-range correlators. {Hi) Elfec- 
tively, higher order contributions of the cluster expansion lead to a screening of long-range tensor 
correlations between two nucleons through the presence of other nucleons within the correlation 
range [24]. For these reasons, we constrain the range of the tensor correlation functions in our 
variational procedure. We use the following integral constraint on the "volume" of the tensor 
correlation functions 

h = ^ drr^ &{r) . (67) 

The constrained energy minimization for the S = 1, T = and the S = l,T = 1 channels with 
dilFerent values of the tensor correlation volume 1^ leads to optimal parameters reported in Table 
121 The optimal parameters for the central correlation functions change only marginally with the 
tensor constraint. Therefore, we adopt a fixed set of parameters for the central correlators given in 
Table E 

The optimal central correlation functions for the Argonne VI 8 potential are shown in Fig. |4l 
In the even channels, the correlation functions decrease rapidly and vanish beyond r « 1.5 fm. The 
central correlators in the odd channels are weaker and of slightly longer range due to the influence 
of the centrifugal barrier. For the tensor correlation functions the constraints on the range are 
important. Fig. [5] shows the triplet-even (a) and triplet-odd (b) tensor correlation functions &(r) 
for different I^ . Because the tensor interaction is significantly weaker for T = 1 than for T = 0, 
the tensor correlator for this channel has a much smaller amplitude. The relevant values for the 
constraint 1^ are therefore smaller for the triplet-odd channel. 

We stress that the range constraint for the tensor correlation functions has an important phys- 
ical and conceptual background. The unitary correlation operator method is used to describe 
state-independent short-range correlations only. Long-range correlations of any kind have to be 
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described by the model space employed in the solution of the many-body problem. By constrain- 
ing the range of the tensor correlators we introduce a separation scale between short-range and 
long-range correlations. The optimal value for tensor constraints cannot be fixed in the two-body 
system alone, but requires input from few-nucleon systems. We will come back to this point in 
SecH 

2.6. Cluster expansion 

The similarity transformation (l23l) of an operator A gives a correlated operator that contains 
irreducible contributions of higher orders in particle number as given by the cluster expansion 

A = C^AC = A"] + JP-'^ + + • • • , (68) 



where A'^"^ denotes the irreducible n-body part 11221] . For a fc-body operator A there will only be 



contributions A^"^ with n > k. For the correlated Hamiltonian we will, therefore, have a one-body 
contribution (from the kinetic energy), a two-body contribution (two-body part of the correlated 
kinetic energy and correlated two-body potential), three-body contributions and so on. 

In practice it will not be possible to evaluate matrix elements of the correlated operators to all 
orders. The importance of the higher-order terms depends on the range of the central and tensor 



correlations 1241 123L 12211 . If the range of the correlation functions is small compared to the mean 
interparticle distance, then three-body and higher-order terms of the cluster expansion are expected 
to be small. In the two-body approximation these higher-order contributions are discarded 

In principle, the higher-order contributions to the cluster expansion can be evaluated systematically 



113 3h . However, for many-body calculations the inclusion of those terms is an extreme challenge 
and we restrict ourselves to the two-body approximation. 

Within the two-body approximation the similarity transformation is still unitary on the two- 
body level, e.g. the eigenvalues of the Hamiltonian are conserved in two-body systems, but it 
is no longer unitary on the many-body level. The energy eigenvalues obtained in exact many- 
body calculations using the correlated interaction in two-body approximation will differ from the 
eigenvalues obtained in exact calculations using the bare interaction. As will be discussed in detail 
in Sec. m we can use exact solutions, e.g., in the No-Core Shell Model framework, to estimate the 
size of the omitted higher-order contributions. 

Technically we can calculate matrix elements of correlated operators in the two-body approx- 
imation of the cluster expansion in any many-body approach using the density matrices p|j|.^^ and 

p*^^^ , , of the uncorrected states 

(0| A |0') ^= ^ p;,;i {k\ At" \m) + pZm A^'^ \mn) . (70) 

km k<l,m<n 

Here the one- and two-body density matrices 

pi;; = ((D|aia™|(D'), p^^^^,, = <0| a^afa^a^ |0') (71) 
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Figure 6: Kinetic (T), potential (V) and total energy (H) obtained with the bare Argonne V18 interaction (left), 
including central correlations (middle) and with central and tensor correlations (right) for the doubly-magic nuclei 
'^He, '^O, and "^''Ca using QHQ. shell-model wave functions. 



are given in a generic single-particle basis \k). Typically we will use harmonic oscillator basis 
states, as the harmonic oscillator basis allows to expand the two-body states in products of relative 
and center-of-mass harmonic oscillator states with the help of the Talmi-Moshinsky transforma- 
tion. 



2. 7. Correlated interaction Vucom 

We define the UCOM interaction Vucom as the two-body part of the correlated Hamiltonian 

H = T + VucoM + V{j'^o^ + ... . (72) 

It contains the contributions of the correlated kinetic energy and of the correlated potential. The 
three-body contribution of the correlated Hamiltonian V^^qj^ has not been evaluated explicitly 
yet. 

If we start from a realistic interaction that is given in an operator representation, e.g. the 
Argonne VI 8 potential, then the UCOM correlated interaction can also be given in operator rep- 
resentation 

Vucom = + O^y^Cr)] , (73) 

p 

where 

Op = {1, ((?r(?2), ql, q'((?i-(?2), h\ L^Cc^j .(^2), (L • S), Si2(f, f), Si2(L,L), ^^^^ 

The dots indicate higher-order contributions of the Baker-Campbell-Hausdorff expansion for the 
tensor transformation that have been omitted. The terms shown above result from a truncation to 
operators of up to fourth order in momentum. For most applications the inclusion of these terms 
is sufficient HH]. 
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In Fig. [6] we show the effect of the correlations on the kinetic and potential energy contributions 
for the Argonne VI 8 interaction. The uncorrelated many -body state is given by the OhD. harmonic 
oscillator configuration (the oscillator parameter is chosen to reproduce the experimental radius). 
Without any short-range correlations the considered nuclei ("^He, '^O, and '^''Ca) are not bound 
at all. Even the potential contributions are repulsive. By including the central correlations the 
potential contributions become attractive but the nuclei are still unbound. Part of the gain by the 
potential contributions has to be paid in form of a larger kinetic energy. Altogether the central 
correlations increase the binding energies by about 20 - 30 MeV per nucleon. With the closed- 
shell trial wave functions used here, contributions from the tensor force are only obtained when 
including tensor correlations. As can be seen this is again a huge effect. The binding energies per 
nucleon increase by about 1 5 - 20 MeV and a total binding of about 4 MeV per nucleon is obtained 
on the OhQ. level using the correlated interaction. 

The existence of an operator representation of Vucom is essential for many-body models that 
are not based on a sim ple oscillator or plane-wave basis. One example is the Fermionic Molecular 
Dynamics model [1351 . l36l . 13711 which uses a non-orthogonal Gaussian basis and does not easily 
allow for a partial wave decomposition of the relative two-body states. Nevertheless, it is possible 
to evaluate the two-body matrix elements of Vucom analytically (radial dependencies are fitted 
by sums of Gaussians), which facilitates efficient computations with this extremely versatile basis 

m. 

As we have emphasized already, the operators of all observables have to be transformed con- 
sistently. The unitary transformation of observables like quadratic radii, densities, momentum 
distributions, or transition matrix elements is straightforward given the toolbox acquired for the 
transformation of the Hamiltonian. The Unitary Correlation Operator Method owes this simplic- 
ity to the explicit state and representation-independent form of the correlation operators. In other 
approaches like the Lee-Suzuki transformation 113 8l . l39l 14011 or the Viowi renormalization group 
method 114 ill , it is not possible to provide a closed form for the effective operators. A discussion of 
effective operators in the Lee-Suzuki approach can be found in [1421. 14311. 

An important feature of Vucom results from the finite range of the correlation functions Ssrir) 
and i?r(r) entering into the generators. Since the correlation functions are of finite range, i.e., 
the correlation operator acts as a unit operator at large r, asymptotic properties of a two-body 
wave function are preserved. This implies that Vucom is by construction phase-shift equivalent 
to the original NN-interaction. The unitary transformation can, therefore, be viewed as a way to 
construct an infinite manifold of realistic potentials, that all give identical phase-shifts. 

It is interesting to observe in which way the unitary transformation changes the operator form 
of the interaction while preserving the phase-shifts. The central correlator reduces the short-range 
repulsion in the local part of the interaction and, at the same time, creates a non-local repulsion 
through the momentum-dependent terms. The tensor correlator removes some strength from the 
local tensor interaction and creates additional central contributions as well as new momentum- 
dependent tensor terms. Hence, the unitary transformation exploits the freedom to redistribute 
strength between local and non-local parts of the potential without changing the phase-shifts. The 
non-local tensor terms establish an interesting connection to the CD Bonn potential, which among 
the realistic potentials is the only one including non-local tensor contributions [|44ll . 
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Figure 7: Two-body density p^"]^^ TUf^iAi) in the S = I, Ms - I, T - Q,Mj - channel of "^He. Left: for the 
uncorrelated trial state |<1)), middle; including central correlations, right: including central and tensor correlations. 
The arrow indicates the orientation of the spin in the two-body channel. 

2.8. Correlated densities 

As we have seen, the inclusion of short-range correlations is essential for obtaining bound 
nuclei when using realistic interactions. However binding energies and spectra provide no direct 
information about these short-range correlations as they are hidden in the correlated or effective 
interactions. 

Other quantities are much better suited to provide insight into correlations. The two-body 
density in coordinate space visualizes the effect of the correlations directly, whereas the one- 
body density in momentum space, the nucleon momentum distribution, allows a comparison with 
experimental data, providing a direct proof for the existence of short-range correlations in nuclei. 

Two-body coordinate space density 

In Fig. |7] we show the two-body density p^]^^, jt^^irn) in the deuteron channel. It has been 
calculated with a OhQ. wave function for "^He. For the uncorrelated wave function the two-body 
density has a maximum at ru = 0, where the potential is strongly repulsive. This defect in the 
wave function is cured by the central correlator which shifts the nucleons apart. For the centrally 
correlated wave function we find now the largest density at |ri2l ~ 1 fm. At this distance the 
potential is most attractive. The tensor force provides attraction if the spins are aligned parallel 
with the distance vector. This is reflected in the two-body density. After applying the tensor 
correlations the density is enhanced at the "poles" and reduced at the "equator". 

The correlated two-body density p^^]^^ jmj^Ai) is calculated from the diagonal part of the two- 
body density in two-body approximation by introducing relative and center-of-mass variables r\2 = 
Xi - X2 and = ^{xi + X2) and coupling the single-particle spins and isospins to total spin S and 
isospin T 

PsMsJMr^^i2,^n) = ^ {^m,^,lmJSMs){^mn,^mt^\TMT)x 

m.,| ,m,j ,m.„ ,m,^ 

'^m\,m,^,m,^,m,^-,m,^,m,i,m,^,m,P^n + 5^12,^12 " ^'^12;^12 + ^^2,^12 " ^'^12) 

(75) 



27 



The dependence on the center-of-mass coordinate is removed by integrating over Xu 

T^sk.TMr^f'n) = Jd'Xn 7^^ljMr(^n,Xn) ■ (76) 

The correlated two-body density is very similar to the results from microscopic calculations 
using the bare Argonne V8 interactions by Suzuki and Horiuchi ll45ll . 



Momentum distributions 

Momentum distributions also directly reflect the existence of short-range correlations. Without 
short-range correlations there would be no high-momentum components in the nuclear many -body 
state. To illustrate this, we calculate the correlated one-body momentum distributions 72(fc) for the 
doubly-magic nuclei "^He and ^^O in two-body approximation. The uncorrected wave functions 
|0) are again OhD. harmonic oscillator configurations, where the oscillator parameters have been 
adjusted to reproduce the experimental radii. The momentum distribution is given by 

nik) = J] m al,,,„,(^)a„,„,„,(^) |0) = J] ('^\ [C^a„\^,,„/^)a„,,,„,(^)C]''' |0) . (77) 

The evaluation of this expression is straightforward but lengthy as it requires an integration over 
single-particle coordinates for correlated wave functions that are expressed in relative and center- 
of-mass coordinates. 

The results for the momentum distributions are shown in Fig. [8l Without short-range corre- 
lations the momentum distributions have no high-momentum components. With only the central 
correlations included, we observe a high-momentum tail, which is almost constant as a function 
of momentum. The contributions from the tensor correlations to the high-momentum tails are re- 
markable. The momentum-distribution from the Fermi surface at about 2 fm"' up to about 4 fm~^ 
is dominated by tensor correlations. Only for very high momenta the central correlations become 
more important. We also observe a strong dependence of the momentum distributions on the range 
of the tensor correlator, especially for smaller momenta close to the Fermi surface. This strong 
dependence on the range of the tensor correlator is caused mainly by the simplified uncorrelated 
wave function used in the present calculation. The eff"ect of long-range correlations induced by 
the long-range part of the tensor force is here only included when using a long-range tensor cor- 
relator. In more realistic calculations these long-range correlations can also be expressed within 
the many-body model space. Other long-range correlations will lead to a further softening of the 
Fermi surface. 

The dominating role of the tensor correlations was also found in other microscopic calculations 



[l46l 1471 . l48n and has been confirmed experimentally by comparing (e, e'pp) with (e,e^n) cross 



exper 

Mm 



sections at high momentum transfer [I491 I50I1 and ip,pp) with (p,ppn) cross sections US 111 . 



2.9. Correlated transition operators 

Short-range correlations have to be considered also for the evaluation of transition matrix ele- 
ments. The transition operators are correlated using the same techniques that have been explained 
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Figure 8: Nucleon momentum distribution T^C^) calculated with closed-shell wave functions for "^He and '*0. Results 

are shown for the uncorrected wave functions ( ), for wave functions including central correlations ( ), 

and including central and tensor correlations with I§ — O.lOfm^^ ( ) and la — 0.20 fm^^ ( ). The gray dots 

indicate the Variational Monte Carlo results from Ref. ll46ll . 



for the Hamiltonian in Sec. 12.41 In case of the radius operator 

^■"^ " A X*-^' ~ ^"""^^ " A2 S 



(78) 



the correlated radius operator in two-body approximation is simply given by 



(79) 



The radius like the quadrupole operator and other long-range operators are only weakly affected 
by the short-range correlations [52]. 

But there are other cases, where short-range correlations are very important. One prominent 
example is the neutrinoless double /3-decay. The effective neutrino potentials appearing here con- 
nect to high momenta and are therefore sensitive to short-range correlations in the many-body 
state. UCOM is nowadays one of the standard approaches for including the short-range correla- 
tions in the transition matrix elements [531. 154 1551. 



2.10. Correlated matrix elements 

For many-body calculations the matrix elements of the correlated interaction Vucom and pos- 
sibly of other operators are needed. For the many-body calculations presented in this review we 
employ the harmonic oscillator basis. Here the two-body matrix elements can be decomposed into 
a relative and a center-of-mass matrix element by means of the Talmi-Moshinsky transformation. 
In the end the relative matrix elements of the correlated interaction 

{n(LS)JMTMT\ Vucom \n'(L' S)JMTMt) = 

<n(L5)7MrM7-|CjCj^HintCf2C, - T^ \n'(L'S)JMTMT) (80) 
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have to be calculated. We use here relative L5 -coupled basis states \n(LS)JM TMj), where n is 
the radial quantum number. The corresponding wave function will be denoted as (pn,L{r), the radial 
wave function as Un^ir), 



{r(LS)JMTMT\n(LS)JMTMT) = (pn,L{r) = 



(81) 



r 



where L is the relative orbital angular momentum, S spin, J total angular momentum, and T 
isospin. The interaction also contains Coulomb and other isospin- and charge- symmetry breaking 
terms so that the matrix elements also depend explicitly on Mj. In the following we will omit the 
M and Mj quantum numbers to simplify the notation. 

In the UCOM approach matrix elements can be calculated in different ways. It is possible to 
expand the correlation operators in the basis states \n{LS)JT) and use the matrix elements of the 
uncorrected interaction. 

An alternative approach is to use the operator representation of Vucom and evaluate the matrix 
elements directly. If one expands the radial dependencies of the individual operator channels in a 
sum of Gaussians, all radial integrals can be calculated analytically. The matrix elements of the 
additional tensor operators contained in Vucom can be given in closed form as well. However, this 
direct approach relies on the truncation of the Baker-Campbell-Hausdorff expansion (1551) . 

This can be avoided by applying the tensor correlations to the basis states which can be done 
exactly. For interactions in operator representation the central correlations will still be applied 
to the Hamiltonian as we have a simple and exact expression for the centrally correlated Hamil- 
tonian (cf. Sec. 12.41) . This approach requires a rewriting of the correlated matrix elements by 
interchanging the order of central and tensor correlation operators using the identity 



The central correlator commutes with Suij, qn) and transforms therefore only ?9(r), see Eq. (|451) . 
The tensor correlator Cn acts on LS -coupled two-body states with L = 7 like the identity operator 
and couples states with L = J +\ with states L = J + I (cf. Sec. 12.31) 



= CjjCj Hint C,-Cn 



(82) 



with the "centrally correlated" tensor correlation operator 



Cn = Cj:CnC, = exp[-ii?(i?+(r)) Si2(?,qn)] ■ 



(83) 



_ (Pn,L{r) 

(KL'S )JT\ Cn \n{LS )JT) = \ cos?Kr)<^„,L(r) 

+ sin ej(r)^n,L(f') 



■,L' 
■,L' 
■,L' 



L = J 



L = J+\ 



7+ 1,L = 7+ 1 



(84) 



where 



0y(r) = 3V/(/+l)W(r)). 



(85) 



Using these relations we can calculate the correlated two-body matrix elements exactly. 



30 



Matrix elements using correlator expansion 

This approach is very general as it works for any unitary transformation of the basis states and 
for any two-body operator. To calculate the matrix element of the correlated interaction 



{n(LS)JMTMT\ Vucom \n'{L' S)JMTMt) = 



{n{LS)JMTMT\ Cjc!, Hint CqC, - Tj^t \n'{L' S)JMTMt) (86) 



we evaluate the correlator in the basis states 

= {niLS)JT\CciC,\n{LS)JT) . 



(87) 



The tensor correlation operator acts as the identity operator in L = 7 channels and couples L = J+l 
with L = J +\ states: 



C 



(JLLST) 



7 dr M*^(r)^ .JWJriunM-ir)) 
/ dr ulj{r) cos 0y(r)^ ^JWJr)Un,L{R-ir)) 
± J dru;-^(r) sin OAr)^ y[Rjr)Un,L{R-{r)) 



L = L = J 

L=L=J±l 

L = J ±l,L = J+l 



(88) 



The correlated matrix elements of the interaction are given with these correlator matrix elements 
as 



{n(LS)JT\yvcoM\n(L'S)JT) 

= {n(LS )JT\ ClCl Hint CnC,- - T^nt \n{L'S )JT) 

= J] C^;!^^'^^\n{LS)JT\n,,An\L'S)JT)C'^^^^^^ ^ ' 

nLfl'L' 

Besides the correlator matrix elements (1881) the matrix elements of the bare interaction are needed 

{n{LS )JT\ Hint \n'{L'S )JT) . (90) 

Here h and h' will now run up to some cut-off Wmax- Convergence is reached only if w^ax is chosen 
large enough. A hard interaction like Argonne VI 8 will connect low momentum states with high 
momentum states up to about 15 fm~^ For typical oscillator constants of lOMeV the summation 
therefore has to extend to n^^x ~ 300. 

Matrix elements for interactions in operator representation 

For a bare potential given in the generic operator representation (l39l) the matrix elements can be 
evaluated using the closed form for the centrally correlated interaction. For the tensor correlations 
only the correlated kinetic energy is given in a closed form without approximations. For the 
potential contributions the tensor correlations will be applied to the wave functions. 

We start with the matrix elements for the local contributions of the form y(r)0 with [r, O] = 
[q^, O] = 0, which includes all operators of the set (l40l) except for the qj. terms. The matrix 
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elements for L = L' = 7 are not affected by the tensor correlations, only the central correlators act 
according to (|54|) . In coordinate representation we obtain 



{n(JS)JT\ClClV(r)OCnCr\n'iJS)JT) = J dru*jir)u„j(r) V(r) {(JS)JT\0\(JS)JT) , (91) 

where V(r) = V{R+{r)) is the transformed radial dependence of the potential. For the diagonal 
matrix elements with L = L' = 7 + 1 we get 

<n(7+l, l)JT\ClCl y(r)0CnC,-|n'(7+l, l)JT) = 



= J dru*jir)Un' 



dr <^i(r) u„Mr) V(r)[ <(7+ 1, l)7r| O |(7+ 1, 1)77) cos%{r) 

+<(7+l, l)7r| 01(7+1, l)7r) sin^ir) (92) 
+<(7+ 1, l)7r| 1(7+1, l)7r) 2 cos Ojir) sin?y(r)] 



with 6j{r) = 9j{R+{r)). Finally, the off-diagonal matrix elements for L = 7 + 1 and L' = 7 + 1 are 
obtained as 



<n(7+l, l)JT\ClCi y(r)0CnC,|n'(7+l, l)JT) = 



drulj,,(r)u„Mr) V(r)[ <(7+l, 1)77] 0|(7+1, 1)JT) cos'?,(r) 

-<(7+l,l)7r|0 1(7+1, l)7r) sin^eXr) ^^^^ 
+<(7+l, l)7r| O 1(7+1, l)7r) cos?,(r) sin'^/r) 
+((7+1,1)771 01(7+1, 1)77) sin?j(r)cos'^j(r)] . 

Apart from the integration involving the radial wave functions, the matrix elements of the operators 
O in LS -coupled angular momentum states are required. Only for the standard tensor operator 
O = Si2(^, p the off-diagonal terms on the right hand side of Eqs. ( 1921 ) and (1931 ) contribute. For all 
other operators in (l40l) the off-diagonal matrix elements vanish, and the above equations simplify 
significantly. 

The effect of the tensor correlator is reflected in the structure of the correlated matrix elements 
(|92l) and (|93l ). It admixes components with AL = +2 to the states. Therefore, the correlated 
matrix element consists of a linear combination of diagonal and off-diagonal matrix elements 
{(LS)JT\ O \(L'S)JT). In this way even simple operators, like or (L ■ §) acquire non- vanishing 
off-diagonal correlated matrix elements (l93l) . 

A closed form is available for the momentum dependent terms of the potential (l39l) . For the 
tensor correlated form of the operator 

V,. = ^[q^V(r) + y(r)q^] (94) 



we obtain 



ClN^rCa = \[qlV{r) + V(r)qj] + y(r)[^'(r)Si2(i^, ^^^^ 
- [qrVm'ir) + i^'(r)y(r)q.]Si2(?,qn) 
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by using Eq. (l58l) . Subsequent inclusion of the central correlations leads to the following expres- 
sion for the diagonal matrix elements with L = U = J in coordinate representation: 



{n(JS )JT\ CjC],V,,CfiC, \n'(JS )JT) = J dr {ii^/r) u„,j(r) 



_ R''(r) 
y(r) W(r) - V'(r)-i^ 

Kirr 



1 V(r) ] 



(96) 



where V'(r) = V'(R-^(r)). As before, the tensor correlator does not affect these matrix elements 
and only the central correlations have to be considered. For the diagonal matrix elements with 
L = L' = J + \ the tensor terms contribute and we obtain 



{n{J+\, \)JT\ ClClW.rC^Cr |n'(7+l, \)JT) 



dr 



V{r) W{r) + V{r) 6' j{rf - V'{r)^^ 



(97) 



with e'j{r) = ej(R+(r)). Likewise, we find 
(niJ+l, 1)771 Cj;C^V,,CnC,|nV±l, 1)^?^) 



± Jdr [<^i(r) ^i(r) - u'*^,ir) u^Mr)] ^^!^''}'^ (98) 



for the off-diagonal matrix elements with L = 7 + 1 and L' = 7 + 1 . 

The matrix elements for the correlated radial and angular kinetic energy can be constructed as 
special cases of the interaction matrix elements discussed above. By setting V(r) = 1 /(IjUrir)) in 
Eqs. (|96l ) to (1981 ) we obtain the matrix elements for the effective mass part of the correlated radial 
kinetic energy The matrix elements of the additional local potential in and the angular 
kinetic energy (ISTI) follow directly from Eqs. (|9TI) to (1931) . 
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3. Similarity Renormalization Group (SRG) 



Unlike the UCOM framework, the concept of the Similarity Renormalization Group (SRG) is 
universal and not tied to the specific correlations relevant in the nuclear many-body problem. The 
renormalization of the Hamiltonian through flow-equations was proposed by Glazek and Wilson 
1 56] in the context of light-front field theory and was further developed by Perry et al. ll57. 581. 
Independently, Wegner [1591. 1601] proposed flow equations for the renormalization of Hamiltonians 
in the context of condensed matter physics. A summary of these developments is given in Refs. 
ll6lL l62n and in Ref.[|63D in this volume. 

Already these initial publications on the SRG contain all formal elements relevant for the 
application in the nuclear physics context, even the specific choice of the generator for the SRG 
flow evolution that will be used in the following was discussed by Szpigel and Perry ll57ll . The first 
application of the SRG for the transformation of a nuclear Hamiltonian was presented by Bogner 
et al. [64j, |65|] and the SRG in connection to the UCOM approach was first discussed in Refs. 



Il66ll67n. 

The general concept of all implementations of the SRG is the transformation of the Hamilto- 
nian to a band- or block-diagonal structure with respect to a specific basis by a continuous unitary 
evolution determined via renormalization-group flow equations. The particular physical system 
and application under consideration determines which basis and generator is used in the flow evo- 
lution. In this respect the SRG approach is very flexible and can be adapted to all kinds of band- 
or block-diagonalizations in any basis of choice [1621. 16811. This flexibility is an advantage of the 
SRG scheme as compared to the UCOM transformation, which is tailored for a very specific type 
of correlations. Moreover, the computational simplicity of the SRG-evolution on the level of ma- 
trix elements opens a clear path towards a consistent evolution of many-body forces beyond the 
level of the two-body cluster approximation (cf. Sec. [2.6[) . Recently, the SRG evolution of a 
nuclear Hamiltonian has been performed on the level of three-body matrix elements %% . thus 
demonstrating the feasibility and power of this scheme. 

3.1. SRG flow equations 



The basic idea of the SRG in the formulation of Wegner [|59l . l6d . |62| . 16111 is to transform the 
initial Hamiltonian H of a many-body system into a diagonal form with respect to a given basis. 
The renormalization group flow equation governing the evolution of the Hamiltonian is given by 



da 



(99) 



where a is a formal flow parameter and Hq, the evolved Hamiltonian. Here we use a general 
operator form of the flow equation in many-body space. The central quantity is the anti-hermitian 
generator rja which determines the physics of the flow evolution. Formally, this is an initial value 
problem with the original Hamiltonian H as initial condition Hq,=o = H. Analogous equations can 
be formulated for the evolution of operators B,), of all observables one is interested in. 



d^a 

da 



(100) 
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Apart from trivial cases, the generator rja will depend on the evolved Hamiltonian Hq, itself. There- 
fore, the flow equation for an observable Bq, cannot be solved independently from the flow equation 
of the Hamiltonian, they have to be solved simultaneously. Formally, we can integrate these flow 
equations defining a unitary operator Uq. of the explicit transformations 

H„ = u;^Hu„ , 

! (101) 

B„ = UlBU„ . 

Note that we chose to define the unitary operator Ua such that the adjoint operator appears on the 
left in the similarity transformation — this is consistent with convention introduced for the UCOM 
transformation (cf. Eq. (1231)). but different from many other discussions of the SRG. From Eqs. 
(IIOII) and (|99l) we obtain a differential equation for the unitary operator U^^, 

d\]„ 

-7^ = -\Jar]a , (102) 

da 

which describes an initial value problem with the trivial initial condition Ua=Q = 1 for the uni- 
tary operator. If the generator would be independent of the flow parameter, e.g. rja = ig with a 
hermitian generator g, then this differential equation could be readily integrated, yielding the stan- 
dard exponential form of the unitary transformation operator U,^ = exp(-iQfg). However, typical 
generators used in the SRG have a non-trivial a-dependence, such that the formal solution for the 
unitary operator does not yield a simple exponential but rather a Dyson series. In practical appli- 
cations it is, therefore, much easier to solve the flow equations (l99l) and (1 1001) directly, without 
reference to the explicit unitary operator. 



The ansatz for the non-hermitian generator rja originally used by Wegner 11591 16011 has a quite 
intuitive structure. It is written as a commutator of the diagonal part of the evolved Hamiltonian, 
diag(HQ,), with the full Hamiltonian Hq,, 

77„ = [diag(HJ,HJ. (103) 

Obviously, the definition of diag(Htt) presumes the choice of a basis — this is the basis with respect 
to which the Hamiltonian shall be diagonalized. By transferring the generator (11031) and the flow 
equation (1991 ) into a matrix representation for this basis, two properties of the flow become evident: 
First, the diagonal form of the Hamiltonian provides a fix point of the flow evolution, since the 
generator rja and thus the right-hand side of the flow equation vanish in this case. Second, the 
off-diagonal matrix elements of the Hamiltonian are continuously suppressed throughout the flow 



evolution, the sum of their squares decreases monotonically 11591 16011. Hence the diagonal form is 
a trivial attractive fixpoint of the SRG flow equation. 

So far this approach is generic and independent of the properties of the particular physical 
system, the Hamiltonian, or the basis under consideration. If considering an A-body system, then 
all the aforementioned relations refer to the operators in A-body space. One of the consequences 
is that even a simple initial Hamiltonian, containing two-body operators at most, acquires up to 
A-body terms in the course of the evolution. For practical applications of the SRG approach in 
the nuclear structure context one, therefore, has to simplify the scheme by confining the evolution 
to two or three-body space, thus discarding higher-order contributions in the evolved interaction. 
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Furthermore, instead of using the diagonal part of the Hamiltonian in the definition of the genera- 
tor, one can use the operator that defines the eigenbasis with respect to which the Hamiltonian shall 
be diagonalized. In this way, we depart from the original goal of Wegner's ansatz to diagonalize a 
Hamiltonian via a flow evolution with a ^ oo and rather aim at the derivation of tamed few-body 
interactions for intermediate values of the flow parameter a that are pre-diagonalized with respect 
to a certain basis. 



A simplified scheme along these lines was suggested by Szpigel and Perry 11571] and applied by 



Bogner and others [|64L 16611 . It confines the evolution to two-body space and uses the generator 

7], = {2nf [Tint, HJ = 2n [f , H„] , (104) 

containing the intrinsic kinetic energy Tint = in the two-body system. The prefactor of 
the commutator is chosen such that the dimension of the flow parameter a is [momentum]""^ or 
[length]"^. It is also common to specify the parameter A = a which has the dimension of mo- 
mentum, instead of the flow parameter a. The square of the two-body relative momentum operator 
can be decomposed into a radial and an angular part, 

f = q^^> ^r = l(4-l + l-4). (105) 

Thus an obvious fix point of the evolution with the generator (11041) is a two-body Hamiltonian Hq, 
that commutes with and L^/r^. Hence, this generator drives the matrix elements of the Hamil- 
tonian towards a band-diagonal structure with respect to relative momentum (q, q') and orbital 
angular momentum (L, L'), i.e., with respect to a partial-wave momentum space representation. 

Though we will only use the generator (11041) in the following, one should note that there are 
many other possible choices for rj^. An evident alternative is to use the single-particle Hamiltonian 
of the harmonic oscillator instead of the kinetic energy in the generator. In this way, the Hamilto- 
nian is driven towards a diagonal form in the harmonic oscillator basis. Using various projection 
operators one can design generators that drive the Hamiltonian towards a block-diagonal struc- 



ture in a given basis 116 8h . This flexibility of the SRG technique holds great potential for further 



refinements and applications of the approach. 

3.2. Evolution of two-body matrix elements 

Starting from an initial two-body Hamiltonian H composed of relative kinetic energy T^i and 
two-body interaction V it is convenient to decompose the SRG-evolved Hamiltonian Hq, in a sim- 
ilar way 

H„ = T,ei + V„. (106) 

All flow-dependence is absorbed in the SRG-evolved two-body interaction Na defined by this 
relation. Rewriting of the flow equation (l99l) using the generator (11041) explicitly for the evolved 
interaction Vq. leads to 



^ = T,ei + VJ = (2//)2 [[T,ei, VJ, T,ei + VJ . (107) 
da 
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Even in this simplified form a direct solution of the operator equation is far from trivial. For 
practical applications we, therefore, work on the level of matrix elements. Though any basis 
in two-body space can be used to define this matrix representation, it is convenient to use the 
eigenbasis of the operator entering into the ansatz for the generator (11041) . In our case, this is the 
4^ operator and it is most convenient to adopt the partial-wave momentum eigenbasis \q{LS)JT), 
where the projection quantum numbers M and Mj have been omitted for brevity. 

In this basis the flow equation (|107l) translates into a set of coupled integro-differential equa- 
tions for the matrix elements 

q') = {q(LS)JT\ V„ \q'(L'S)JT) . (108) 
In a generic form, the resulting evolution equation reads: 

^Va{q,q') = -{q^ -q'^f Va(q,q') 

^ (109) 
+ 2^1 dQQ^ {q^ + q'^ - 2Q^)VMQ)ya{Q,q') . 

For non-coupled partial waves with L = U = J, the matrix elements entering this equation are 
simply 

VM,q') = yf^\q,q')- (110) 

For coupled partial waves with L,L' = 7+1, the Va{q,q') are understood as 2 x 2 matrices of 
the matrix elements for the different combinations of the orbital angular momenta L = J - I and 

Each non-coupled partial wave and each set of coupled partial waves evolves independently of the 
other channels of the interaction. This is a direct consequence of the choice of the generator — the 
evolution towards a diagonal in momentum space is done in an optimal way for each individual 
partial wave. 

As mentioned earlier, analogous evolution equations have to be solved for all observables in 
order to arrive at a consistent set of effective operators. The evolution of these operators, e.g. the 
multipole operators necessary for the evaluation of transition strengths or the one-body density 
operators employed for the computation of the momentum distribution, is coupled to the evolu- 
tion of the Hamiltonian via the generator rja. Hence we have to solve these evolution equations 
simultaneously. 

An alternative approach is to determine the matrix elements of the unitary operator Uy, explic- 
itly by solving (11021) . The evolved matrix elements of all observables can then be obtained by 
a simple matrix transformation using the same unitary transformation matrix. In the case of the 
momentum-space partial-wave matrix elements of the unitary transformation operator, 

Ui''^'^''^\q, q') = {q{LS)JT\ U„ \q'{L'S)JT) , (1 12) 
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the operator equation (|102|) leads to a coupled set of integro-differential equations 



1 r> 

— UM^q') = llx dQQHq'^ - Q^) U^(q,Q)Va(Q,q') , (113) 
da J 

where we assume that the evolution equation (11091 ) is solved simultaneously providing the V^iq, q'). 
The generic notation defined in (II 101) and (II 1 II) for non-coupled and coupled partial waves, respec- 
tively, applies here as well. This differential equation provides direct access to the matrix elements 
of the unitary operator, which maps the initial operators onto any particular point of the flow 
trajectory. 

3.3. Evolved interactions and wave functions 

The concept of the SRG transformation becomes very transparent when looking at the flow 
evolution of the momentum- space matrix elements of the SRG-transformed interaction V^. In 
Fig. |9]we show the matrix elements Va^^ ^^\q, q') obtained for the Argonne V18 potential in the 
three most important partial waves: the partial wave, i.e. matrix elements Va^^^^\q,q'), the 
1 partial wave, i.e. yr°'°^(^, q'), and the i - ^Dj partial wave, i.e. V[^^'^^^\q, q'). We start 
with the matrix elements of the initial Argonne VIS potential at a = Ofm"* and display snapshots 
of the SRG evolution at a = 0.001 fm\ a = 0.01 fm"^, and a = 0.04 fm^ 

The initial matrix elements show the characteristic features that are responsible for the emer- 
gence of strong correlations in the many-body system: The strong off-diagonal matrix elements 
that couple low-momentum components with high-momentum components of the wave function. 
In the and i partial waves these off-diagonal high-momentum matrix elements are generated 
by the short-range repulsion of the Argonne VI 8 potential, while the off-diagonal matrix elements 
in the i - ^Di partial wave are solely due to the tensor interaction. 

Already in the early phase of the flow evolution, i.e. for a < 0.01 fm"^, the matrix elements 
far off the diagonal in the S -wave channels are suppressed quickly. The plateau of positive high- 
momentum matrix elements is pushed towards the zero-plane and the negative low-momemtum 
matrix elements are enhanced. Later in the flow evolution, the residual high-momentum matrix 
elements are pushed towards the diagonal and the matrix elements in the low-momentum are 
further enhanced. For the tensor-dominated \ - ^Di partial wave the matrix elements far-off 
the diagonal are depleted successively. Over all, the SRG-evolution with the generator (11041) 
leads to a transformed interaction with band-diagonal matrix elements in momentum space. The 
pre-diagonalization of the interaction and thus the Hamiltonian matrix elements through the SRG 
transformation is evident. 

In order to assess the effect of the flow-evolution on the structure of wave functions, we again 
consider the deuteron ground state as an example. Using the SRG-evolved momentum-space 
matrix elements we solve the two-body problem in the deuteron channel numerically. The de- 
pendence of the resulting coordinate- space wave functions on the flow parameter a provides an 
intuitive picture of the effect of the SRG evolution in coordinate space. This is illustrated in Fig.fTOl 

The deuteron wave function for the initial interaction shows the signatures of strong short- 
range central and tensor correlations, i.e., the suppression of the wave function at small inter- 
particle distances and the presence of the D-wave admixture, respectively. Already early in the 



38 



^^0 ^Si ^Si-^Di 




Figure 9: Snapshots of the momentum-space matrix elements (in units of MeV fm^) of the SRG-evolved Argonne 
V18 potential (charge independent parts only) in the '5o, i, and \-^Di partial waves (from left to right) for the 
flow parameters a - 0.0 fm**, 0.001 fm"*, 0.01 fm"*, and 0.04 fm'* (from top to bottom). 
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Figure 10: Snapshots of the deuteron wave function obtained from the SRG-evolved Argonne V18 potential for the 
flow parameters a - 0.0 fm^, 0.001 fm'*, 0.01 fm'*, and 0.04 fm'* (from top to bottom). The main panels show the radial 

wave functions (piir) for L = ( ) and L - 2 { ). The 3D plots show the corresponding momentum-space 

matrix elements in the i partial wave for orientation (cf. Fig. |9]l. 
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flow evolution, i.e. for a < 0.01 fm'*, the short-range dip in the 5 -wave component is removed — 
this effect is connected to the suppression of the plateau of high-momentum matrix elements in 
the ^Si partial wave. As a result the deuteron wave function obtained with the SRG-evolved 
Argonne VI 8 potential for a = 0.01 fm"^ has lost any signature of a strong short-range repulsion 
in the interaction. The total strength of the D-wave admixture is reduced and pushed towards 
larger inter-particle distances. Thus, the short-range wave function evolves from being dominated 
by short-range central and tensor correlations to an almost uncorrelated pure S -wave. The long- 
range behavior, all asymptotic properties, and the deuteron binding energy are not affected by the 
SRG-evolution. 

The behavior of the matrix elements and of the two-body wave functions highlights the relation 
between the UCOM transformation and the SRG evolution regarding the pre-diagonalization of 
the Hamiltonian and the description of short-range correlations. Both approaches describe the 
same physics. 

3.4. UCOM from an SRG perspective 

Because the UCOM and the SRG transformations have the same effect on the matrix elements 
and wave functions, one might ask for the connection of both approaches on the underlying for- 
mal level. The properties of both unitary transformations are governed by their generators: the 
dynamical generator -iria in the case of the SRG and the static generators g, and gn in the case of 
the UCOM transformation. 

There is a non-trivial relation between these generators that provides an insight into the formal 



relation and the differences of the two approaches 116611 . This becomes evident by evaluating the 
SRG generator (11041) for a typical nuclear Hamiltonian in two-body space. We assume a simplified 
local two-nucleon interaction composed of a central, a spin-orbit and a tensor part. The operator 
for this interaction is given by 

Y = Y,VpiT)0, (114) 

p 

with Op e {1, ((Ti •(T2), (L • S), Si2(;, ;)} <S> {1, (ti -72)}. By evaluating the commutator defining the 
SRG generator (I104|) explicitly for a = using this interaction operator we obtain 



^{qr S(r) + S(r) q,) + 0(r) Si2(?, 4ci) 



(115) 



The operator- valued functions S(r) and 0(r) contain the radial dependencies of the different terms 
of the interaction 

If the functions S(r) and 0(r) were functions of the relative distance r alone, then the structure 
of the initial SRG generator irjo would be identical to the UCOM generators g,. and g£i that were 
constructed based on the physical picture of central and tensor correlations [;66i1 . The explicit 
dependence of S(r) on the operator set O^^ indicates that the SRG transformation acts differently 
in different partial waves. In the UCOM terminology, this dependence encodes a partial- wave 
dependence of the central correlation functions. 
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This formal connection shows that both approaches address the same physics of short-range 
correlations, although starting from quite different backgrounds. Moreover, it proves that the set 
of UCOM generators covers the most relevant terms. Although there are other operators appearing 
in the initial interaction, e.g. the spin-orbit operator, they do not require separate generators — their 
effect on the correlations is absorbed in the operator-valued function S(r). 

Regarding the partial-wave dependence of the correlations, the standard formulation of the 
UCOM approach uses a simplified picture. The correlation functions s(r) and 7?(r) are chosen to 
depend on spin S and isospin T only, they do not depend on orbital and total angular momentum. 
Formally, one could drop this restriction and work with separate correlation functions for each 
partial wave and thus mimic the flexibility of the SRG generator. 

A more fundamental difference between the UCOM and the SRG transformations results from 
the fact that SRG uses a dynamical generator, whose operator structure changes throughout the 
flow evolution, whereas UCOM is based on a static generator. Though the UCOM generator and 
the SRG generator share the same basic operator structure at a = 0, the SRG generator acquires a 
more complicated form involving higher-order momentum and momentum-dependent tensor op- 
erators at later stages of the evolution. Therefore, the SRG generator is more flexible and adapts 
to the behavior of the matrix elements during the flow evolution, leading to a non-trivial flow 
trajectory in an operator-space representing the generator. The UCOM transformation, in con- 
trast, consists of a one-step transformation along a linear trajectory confined to a subspace of the 
operator- space spanned by the SRG generator. It is, therefore, not surprising that the matrix ele- 
ments of the UCOM-transformed interaction do not exhibit the same perfect band-diagonal struc- 
ture as the SRG-evolved interaction. However, the leading operator contributions to the generator 
are also present in the UCOM approach and allow for an efficient pre-diagonalization. 

The dynamic nature of the SRG generator is also the reason, why the optimal UCOM correla- 
tion function cannot be determined directly from Eqs. (II 161) . We have to consider the whole SRG 
flow trajectory up to a certain value of a to extract meaningful UCOM correlation functions. One 
option to do so is discussed in the following section. 

3.5. UCOM correlation functions extracted from SRG 

As an alternative to the variational determination of the UCOM correlation functions R+ir) 
and ??(r) discussed in Sec. 12.51 we can use the SRG approach to generate optimized UCOM cor- 
relation functions. Our aim is to construct a UCOM transformation that uses the result of the 
SRG-evolution of a Hamiltonian in two-body space from or = to a fixed finite value of the flow 
parameter to determine R+ir) and ??(r). In contrast to the dynamical SRG evolution, the UCOM 
correlations have to map the initial Hamiltonian onto the evolved Hamiltonian for a specific a 
using a single explicit unitary transformation. Obviously, the correlation operator C is not flexi- 
ble enough to allow for an exact mapping of all matrix elements in all partial waves — even if we 
would allow for different correlation operators C for each partial wave. 

One could consider an approximate mapping of the matrix elements as one possible scheme 
to determine the UCOM correlation functions. Here, we use a different strategy, which is rooted 
in the interpretation of the UCOM transformation as a tool to imprint short-range correlations into 
the many-body state. Instead of considering the initial and evolved two-body matrix elements, we 
consider two-body eigenstates of the initial and evolved Hamiltonian for different partial waves. 
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The optimal UCOM correlation functions are then required to map a selected two-body eigenstate 
of the SRG-evolved Hamiltonian onto the corresponding eigenstate of the initial Hamiltonian. 
This wave-function mapping defines the so-called SRG-generated UCOM correlation functions. 

The procedure for the construction of SRG-generated UCOM correlation functions consists 
of three steps: (?) We solve the SRG evolution equations for a given initial interaction up to a 
flow parameter a, obtaining the momentum space matrix elements Va(q, q') for a certain partial 
wave, (zz) Using the evolved matrix elements the two-body problem is solved, leading to a set 
of coordinate-space wave functions, (zzz) The UCOM correlation functions R+{r) and §{r) are 
determined such that they map a selected two-body eigenstate of the SRG evolved interaction onto 
the corresponding two-body state of the initial interaction in the respective partial wave. 

The steps (z) and (zz), i.e., the evolved momentum- space matrix elements and the wave func- 
tions of the corresponding two-body eigenstates, respectively, have already been illustrated for the 
deuteron channel in Sec. 13.31 Step (zzz) is discussed in the following. 

Consider two eigenstates 1^*"-') and 1^''^"'') with the same energy eigenvalue resulting from the 
solution of the two-body problem for the initial and the SRG-evolved potential, respectively, in a 
given coupled or non-coupled partial wave. We can define a UCOM correlation operator C that 
maps the two states onto each other 

|<^(0)) = C|(^('^>) = CnC,|^('^)). (117) 

Based on this formal definition we can derive equations that determine the correlation functions 
R-{r) and d-{r) that characterize the correlation operator. 

For non-coupled partial waves with L = J only the central correlator is relevant. With the 
two-body solutions 

|.^(°)) = i0(°)(LS)yr) 

(118) 

l^^''') = \(f>^"\LS)JT) 

for the initial and the SRG-evolved interaction, respectively, we obtain from (II 171) and (l32l) a rela- 
tion connecting the known radial wave functions 0^"^(r) and (p^^^r) via a yet unknown correlation 
function R-{r): 

0(O)(^) = (p^"\R-(r)) . (119) 

r 

Here and in the following we assume real- valued wave functions. The relation (II 191 ) can be viewed 
as a differential equation for the correlation function R-(r). After formal integration we arrive at 
an implicit integral equation for R-{r) 

[R_(r)f = 3 d^f .Z.r: Lo ' (120) 



which can be solved easily in an iterative fashion. We end up with a discretized representation of 
the correlation function R-{r) for the partial wave under consideration. By construction it maps a 
selected SRG-evolved two-body state onto the corresponding initial state. In general, R-{r) will 
depend on the pair of states, e.g. the ground states or a pair of excited states, we have selected. We 
will show later on that this dependence is very weak. 
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For coupled partial waves with L = J - 1 and L' = J + I central and tensor correlators act 
simultaneously. Using the two-body eigenstates 



1^^°)) = \cf>f(LS)JT) + \cf>f,\L'S)JT) 
\^^"^)= \cpf{LS)JT)+\cpf{L'S)JT) 

of the initial interaction and the evolved interaction, respectively, we can extract a unique set 
of central and tensor correlation functions. After multiplying the mapping equation (II 171) with 
{r{LS)JT\ and {r{L'S)JT\ , respectively, and using Eq. (iJTl) . we obtain a system of coupled equa- 
tions 

0f (r)\ _ R-{r) r-— I cos OAr) sin ej{r)\ lcpf{R.{r))\ 

~ ~^ ^ ^ \- sin ej{r) cos OAr)) V^{R.{r))) ' ^^^^^ 

from which the correlation functions R-{r) and §{r) can be determined. 

Because the central correlation function acts on both orbital components in the same way 
and because the transformation matrix in (11221) has to be unitary, we can determine the central 
correlation function R-{r) independently of the tensor correlations function ??(r). By considering 
the sum of the squares of the two orbital components we obtain from (11221) the identity 



VR-ir)f 

r2 



[4 W] + [-^L' ('■)] = —^R'-^r) ([4 (^-(r))]' + [4?(i?-(r))]0 . (123) 



which corresponds to (II 191 ) for the non-coupled case. The correlation function R-{r) can then be 
determined iteratively from the integral equation 



{R-{r)r = 3\ > : 7; . (124) 

Once R-{r) is known, the system (11221) reduces to a set of two nonlinear equations for 9j{r) = 
3 V^TO/TT) &(r), which can be solved numerically for each r. 

In practice this mapping scheme can be easily implemented using discretized wave functions. 
Typically, the SRG evolution of the Hamiltonian in two-body space for a given partial wave is per- 
formed on a sufficiently large grid in momentum space. Using the discretized momentum- space 
matrix elements of the evolved Hamiltonian we solve the two-body problem for the respective 
partial wave on the same momentum- space grid. The ground state wave-functions are then trans- 
formed to coordinate space, where the mapping equation (I122|) is solved. In this way, we obtain 
discretized correlation functions R-{r) and, by numerical inversion, R+(r) as well as i?(r) for each 
partial wave. In contrast to the UCOM correlation functions determined variationally, there are no 
parameterizations of the correlation functions necessary, which might induce artifacts due to their 
limited flexibility. 

3.6. SRG-generated UCOM correlation functions 

As an example for the determination of UCOM correlations functions through a mapping of 
SRG-evolved wave functions, we again consider the Argonne VI 8 potential. Note, however, that 
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the mapping procedure is completely generic and can be used with any other interaction, be it local 
or non-local. 

In order to stay within the framework set by the variational UCOM correlators, we do not 
consider separate correlation functions for each partial wave, which could be done easily, but only 
distinguish different channels of two-body spin S and isospin T. As in the variational scheme, we 
optimize the correlators for a given (5, r)-channel using the lowest 7 partial wave, since the low-7 
and thus low-L partial waves are affected most by short-range correlations. For partial waves with 
higher J and L the impact of short-range correlations and of the UCOM transformation is reduced 
due to the angular momentum barrier, thus the non-optimal correlators for these partial waves do 
not have a big impact. 

For determining the central correlation functions R+(r) in the spin-singlet channels, we use 
the partial wave for the singlet-even (5 = 0, T = I) and the ^Pi partial wave for the single- 
odd (S = 0, r = 0) channel. The central and tensor correlation functions in the triplet-even 
(5 = 1, r = 0) channel are extracted from the deuteron solution in the coupled i - ^Di partial 
wave. For the triplet-odd channel (5 = I, T = 1) we encounter the same ambiguity as in the 
variational treatment: the lowest possible orbital angular momentum allowed by antisymmetry is 
L = 1 for which J can be 0, 1, or 2 and only for 7 = 2 the tensor correlator does contribute. One 
possible recipe for handling this channel is to use only the coupled ^Po - partial wave to fix the 
triplet-odd central and tensor correlation functions, as done in Ref. 11671] . The central correlation 
functions obtained in this way are not well adapted for the lower-7 partial waves. Therefore, a 
scheme that includes all possible 7 for the determination of the central correlation function seems 
more appropriate. Thus, we determine the central correlation function through a mapping of wave 
functions for a pseudo interaction obtained by averaging the ^Pq, ^Pi, and ^P2 partial waves with 
a relative weight 27 -i- 1 . This recipe comes closest to the energy-average used for the variational 



determination in this channel (cf. Sec. 12.51) . 

In each of the partial waves we use the energetically lowest pair of states, which is bound in 
the case of the triplet-even channel and unbound otherwise, to determine the correlation functions 
via the mapping (II 171) . In principle one could use any other pair of two-body states obtained for 
the initial and the SRG-evolved Hamiltonian with the same energy. As was shown in Ref. 11671] 
the correlation functions do not change significantly when using one of the low-lying excited two- 
body states instead of the ground state. 

A crucial advantage of the SRG-generated correlation functions is that there is no need for 
artificial constraints to control the range of the correlations functions — the SRG flow-parameter a, 
which enters through the evolved two-body eigenstates, is the only control parameter. In contrast to 
the ad-hoc integral constraints formulated for central and tensor correlations functions in Sec. 12.51 
the flow parameter a is a physically motivated control parameter that enters the central and tensor 
correlation function in a consistent way. 

The dependence of the correlation functions on the flow-parameter a is illustrated in Fig. \TT\ 
for the central correlation functions and in Fig. [12] for the tensor correlation functions obtained 
for the Argonne VI 8 potential through the SRG mapping. Evidently, the over- all range of the 
correlation functions is directly controlled by the flow-parameter a: Larger a result in correla- 
tion functions with longer range. This is in-line with our observations on the evolution of the 
momentum- space matrix elements and the two-body wave functions. Initially, the SRG flow- 
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Figure 1 1 : Illustration of the dependence of the SRG-generated central correlation functions R+(r) - r for the Argonne 
V18 potential on the flow parameter a for the different spin S and isospin T channels. The curves correspond to 
a = 0.03fm'* ( ), a = 0.04fm'* ( ), a = 0.06fm'* ( -), and a = O.OSfm'' ( ). 




Figure 12: Illustration of the dependence of the SRG-generated tensor correlation functions ??(r) for the Argonne 
V18 potential on the flow parameter a for the different spin S and isospin T channels. The curves correspond to 
a = 0.03fm'* ( ), a = 0.04fm'* ( ), a = 0.06fm'* ( -), and a = O.OSfm'' ( ). 
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evolution affects only the high-momentum matrix elements and thus the short-distance behavior 
of the wave functions. Throughout the flow evolution, i.e., with increasing a, the wave functions 
are modified at increasingly larger distances, which results in an increasing range of the associated 
correlation functions. This localized action of the SRG transformation on coordiante-space wave 
functions is also responsible for the fact that the SRG-generated correlation functions automati- 
cally have finite range. This property is not imposed by the mapping scheme, but results from the 
structure of the two-body wave functions alone 3. 

Closer inspection of the structure of the correlation functions in Figs. [TT] and [12] reveals an 
interesting new aspect as compared to the correlation functions discussed in Sec. [21 In the even 
channels the central correlation functions R+(r) - r exhibit a sign-change at r « 1.1 fm. At shorter 
ranges R+(r) - r is positive, indicating an outward shift in a transformed two-body wave func- 
tion, and turns negative at larger distances, inducing an inward shift in a transformed wave func- 
tion. Pictorially speaking, the UCOM transformation attempts to exploit the attractive parts of 
the central potential in the even channels by moving probability amplitudes from small and large 
inter-particle distances into the attractive region. The triplet-even tensor correlation function &(r) 
exhibits a similar structure, though the negative contribution is much weaker than the positive part. 
In the odd channels the correlation functions do not show this sign change, which can be explained 
by the lack of a sufficiently strong attraction in the central interaction. 

The details of the ^-dependence are different for the different types of correlation functions. 
For the dominant central correlators in the even channels, the short-range positive component is 
practically independent of a in the range covered in Fig. [TT] — these are the generic short-range 
correlations induced by the strong short-range repulsion of the Argonne VI 8 potential. Only the 
negative long-range part shows a sizable a-dependence affecting its range and strength. All other 
correlation functions show a smooth increase of the over-all range with increasing a, only a very 
short distances the curves are independent of a. 

The direct comparison of the SRG-generated correlation functions with the correlators deter- 
mined variationally in Sec. I2.5l is also quite instructive. In Figs. [T3and[T3|we compare the central 
and tensor correlations functions, respectively, of the SRG-generated correlator for a = 0.04 fm'* 
and the standard variational correlator with constraints = 0.09 fm^ for the tensor correlator in 
r = 0, and Iff = -0.03 fm^ for the tensor correlator inT = I. Both sets of correlators yield approx- 
imately the same ground-state energy of "^He in a No-Core Shell Model calculation, as discussed 
in Sec.[S 

In the even channels, the short-range parts of the central correlation functions of both sets 
agree very well — another indication that these dominant short-range correlations are truly generic 
and independent of the methodology used to determine the correlation functions. The negative 
contributions in R+(r) - r appearing in the SRG-generated correlators is absent in the variational 
correlators, simply because the parameterizations used for the latter did not allow for such a struc- 
ture. The triplet-odd central correlators also agree very well, thus providing additional justification 
for the treatment of this channel in the mapping procedure. 



'We note that the short range of the correlation functions depends on the initial interaction. If the interaction 
is such that the SRG evolution affects also long-range components of the wave function, then the UCOM correlators 
obtained by the mapping will be long-ranged as well. This is the case for the chiral N3LO interaction fl, for example. 
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Figure 13: Comparison of the central correlation functions R+{r) - r for the Argonne V18 potential obtained via 

variational optimization with constraints /^'"' = 0.09 fm"^, and = -0.03 fm"' ( ) and via SRG-mapping with 

a = 0.04fm4 ( ). 
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Figure 14: Comparison of the tensor correlation functions »?(r) for the Argonne V18 potential obtained via variational 

optimization with constraints tf^^ - 0.09 fm^, and /^''' = -0.03 fm^* ( ) and via SRG-mapping with a - 

0.04 fm"* ( ). 
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In the singlet-odd channel the central correlation functions exhibit the same shape, but the 
variational correlator is suppressed due to the explicit range constraint. However, the variational 
correlation function agrees very well with the SRG-generated central correlators for smaller values 
of a. The tensor correlations functions determined in the variational scheme are also subject to 
explicit range constraints that affect their shape. Over-all the correlation functions &(r) of both 
sets are similar, but the agreement is not as good as for the channels without ad hoc constraints. 

For the following discussion, we will identify the SRG-generated UCOM correlation func- 
tions and the resulting UCOM-transformed potential with the abbreviation "UCOM(SRG)". The 
UCOM correlation functions determined from a variational calculation (cf. Sec. 12.51) are termed 
"UCOM(var.)". Finally, the purely SRG-transformed interactions are labelled "SRG". 

3.7. Comparison of matrix elements 

To conclude the discussion of the transformed interactions resulting in the UCOM and in the 
SRG scheme, we consider the momentum-space matrix elements once more. In Fig. [15] we com- 
pare the matrix elements in the dominant 5 -wave channels obtained for the bare Argonne VI 8 
potential, the UCOM-transformed Argonne VI 8 potential using correlation functions obtained 
variationally as well as via the SRG mapping, and the pure SRG-evolved interaction. 

All similarity transformed interactions show a strong suppression of the off-diagonal matrix 
elements, i.e., a decoupling of low-momentum and high-momentum states, and an enhancement 
of the low-momentum matrix elements. This leads to a significant improvement of the convergence 
in No-Core Shell Model calculations for light nuclei, as will be discussed in detail in Sec. HI There 
are, however, distinctive differences in the behavior of the matrix elements for the and the \ 
partial waves in the high-momentum sector. 

The SRG evolution leads to a transformed interaction with a perfect band-diagonal structure in 
momentum space, i.e., in the high-momentum regime the matrix elements drop to zero rapidly with 
increasing distance from the diagonal. For the UCOM-transformed interaction, the domain of non- 
vanishing high-momentum matrix elements extends further out. In the case of the UCOM(var.) 
correlation functions there is a plateau of non-vanishing matrix elements in the high-momentum 
sector, which falls off slowly when leaving the diagonal. As a result, the band-diagonal structure 
is far less pronounced than in the case of the SRG-evolved interaction. For the UCOM(SRG) 
interaction there appears a broad band of non-vanishing high-momentum matrix elements. The 
far-off diagonal matrix elements outside of this band are more suppressed than for the UCOM 
with variationally determined correlators, but compared to the SRG-transformed interactions the 
band is significantly broader. 

The apparent differences between the SRG- and the UCOM(SRG)-transformed matrix ele- 
ments show that, despite the construction of the UCOM correlation functions using input from 
the SRG evolution, the transformed interactions are very similar for ^q^ + q'^ < 2fm"\ but are 
quite different above. The origin of this differences is the limited flexibility of the UCOM gen- 
erator, which is determined from the eigenstates with the lowest energy and does not allow for a 
perfect band-diagonalization of the high-momentum matrix elements. However, it does allow for 
a decoupling of low- and high-momentum modes, which will be important for the convergence of 
shell-model-like many-body calculation as discussed in the following section. 
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Figure 15: Momentum-space matrix elements (in units of MeV fm^*) of the UCOM and SRG-transformed Argonne 
V18 potential in the '^o, ^^i, and ^^i-^Di partial waves. First row: initial Argonne V18 matrix elements (charge 
independent terms only). Second row: UCOM transformed matrix elements using UCOM(var.) correlation functions 
for 7^'°^ = 0.09 fml Third row: UCOM transformed matrix elements using UCOM(SRG) correlation functions for 
a - 0.04 fm^ Fourth row: SRG-transformed matrix elements for a - 0.03 fm^. 
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4. No-Core Shell Model 



The No-Core Shell Model (NCSM) is a powerful and well established many-body technique 
that has been used successfully in a wide range of nuclear structure calculations in light nuclei 



11701 17 11.17211. It provides a perfect framework for assessing the properties of nuclear interactions, 
both from the technical perspective, e.g., regarding their convergence properties, and from the 
experimental view, e.g., regarding the agreement of predicted observables with experiment. 

In this section we discuss NCSM calculations using the UCOM interaction with correlation 
functions determined variationally (cf. Sec. [21) and the UCOM interaction using SRG-generated 
correlation functions (cf. Sec.O). For comparison we also show results with the SRG interaction. 
All UCOM and SRG interactions are derived from the Argonne VI 8 interaction. The matrix 
elements are calculated without approximation as explained in Sec. 12.101 and Sec. l3.2l including all 
electromagnetic and charge dependent terms. 

Using exact calculations in the three- and four-body system, we investigate the convergence 
properties and the role of induced three-body interactions of the UCOM and SRG interactions as 
a function of tensor correlation range or flow parameter. The binding energies in the three- and 
four-body system map out the so-called Tjon-line. By choosing particular values for the tensor 
correlation range and the flow parameters, respectively, we obtain UCOM and SRG two-body 
interactions that provide binding energies very close to the experimental results in the three- and 
four-body system. 

The UCOM and SRG interactions selected by this choice are then used in NCSM calculations 
of ^He, ^Li and ^Li. Although we are not able to reach full convergence, binding energies can 
be obtained by extrapolation. A better convergence behavior is found for the excitation energies. 
The spectra provide some insight about the spin dependence of the UCOM and SRG interactions. 
Additional hints are provided by electromagnetic properties like radii, magnetic dipole moments 
and electric quadrupole moments. 

4.1. Benchmarking Vucom in ab initio few -body calculations 

For the NCSM calculations in the three- and four-body system we use the ManyEff code by 



Petr Navratil jllll . It employs a translationally invariant oscillator basis in Jacobi coordinates. The 
model space is defined by the oscillator frequency hO. and the total number A^^ax of excitations 
with respect to the O^Q configuration. The interaction is provided in form of relative harmonic- 
oscillator matrix elements. In the standard NCSM approach an effective interaction adapted to the 
model space (N^ax, is derived using the Lee-Suzuki transformation. For the calculations pre- 
sented here we directly use the "bare" UCOM and SRG matrix elements, which do not depend on 
the model space size N^nax- Within this procedure, the NCSM provides a variational approach for 
the energy. The energy eigenvalues will converge from above to the exact solution for sufficiently 
large model spaces. The converged results should also be independent of the oscillator frequency 
h£l, although the rate of convergence will be different for different oscillator frequencies. 

Before presenting results with UCOM and SRG interactions we illustrate that it is not possi- 
ble to reach full convergence with the bare Argonne VI 8 interaction, even within the huge model 
spaces possible for and "^He. As can be seen in Fig.[T6]it is already very hard to obtain a bound 
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Figure 16: NCSM calculations with the bare Argonne V18 interaction for and "^He. The ground state energy is 
calculated in model space with A^max = 0, 2, . . . , 40 for and with A^max = 0, 2, . . . , 18 for "^He. The exact binding 
energies jl^l are indicated by horizontal Unes. 
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Figure 17: NCSM calculations for the ground-state energy of H in model spaces with A^max — 0,2,.. .,40. The 
tensor correlation range for the UCOM(var) interaction is = 0.09 fm^. The flow parameters are a - 0.04 fm'* 
for the UCOM(SRG) and a — 0.03 fm"* for the SRG interaction. The experimental binding energy is indicated by a 
horizontal line. 



nucleus. One can also observe that the lowest energies are obtained for very large oscillator fre- 
quencies. With these narrow oscillator wave functions it becomes eventually possible to explicitly 
describe the short-range correlations for very large A^max- 

We can compare these results with calculations with the UCOM and SRG interactions as shown 
in Fig.[l7]and Fig.[T8l We use our "standard" choices given in (11251) for the tensor correlation range 
and flow parameters a as explained later in this section. For all interactions we find a bound 
minimum already in the OhO. space. The minima are at oscillator frequencies that correspond 
roughly to the experimental sizes of the nuclei. With increasing model-space size N^i,^ we observe 
a fast convergence for all interactions. In direct comparison, the UCOM(SRG) and SRG inter- 
actions converge faster than the UCOM(var.) interaction with respect to both, model space size 
A^max and oscillator frequency hO.. It is important to note that the converged energy, while being 
close to the experimental binding energy, is lower than the exact result for the bare Argonne VI 8 
interaction by about 0.7 MeV in case of and by about 4.1 MeV in case of "^He. This overbind- 
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Figure 18: NCSM calculations for the ground-state energy of ''He in model spaces with A^max - 0, 2, . . . , 18. The 
tensor correlation range for the UCOM(var) interaction is = 0.09 fm^. The flow parameters are a - 0.04 fm'^ 
for the UCOM(SRG) and a - 0.03 fm"* for the SRG interaction. The experimental binding energy is indicated by a 
horizontal line. 
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Figure 19: NCSM calculations for the ground state of "'He with UCOM(SRG) interactions at different flow parameters 
(a - 0.01, 0.04, 0.08 fm"*) in model spaces with A^max = 0, 2, . . . , 18. The horizontal line indicates the experimental 
binding energy. 



ing of UCOM and SRG interactions with respect to the bare interaction is caused by the missing 
three- and four-body contributions - UCOM and SRG interactions are only calculated in two-body 
approximation. This point will be discussed in detail in Sec. 14.21 

The convergence pattern and the converged energy depends on the parameters of the tensor cor- 
relation range or the flow parameter respectively. This is illustrated in Fig.[T9]for the UCOM(SRG) 
interaction using three diff"erent flow parameters. With increasing flow parameter the calculations 
converge faster and to a lower energy. This is analyzed in detail in Fig. [201 where we compare 
the energy minima in the OhQ model spaces and the converged energies as a function of 1^ or a. 
The OhQ results are getting closer to the converged results with larger correlation ranges or flow 
parameters, i.e., the interactions become "softer". Note that the converged UCOM(var.) results 
decrease monotonically with increasing tensor correlation range, whereas the UCOM(SRG) and 
SRG interactions both show a minimum in the converged energy for flow parameters or « 0. 10 fm"^. 
This indicates that the tensor correlation range and the flow parameter in UCOM(SRG) and SRG 
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Figure 20: "'He ground-state energy calculated in QhQ. space (dotted lines) and converged NCSM results (solid lines) 
as a function of correlation range or flow parameter. The exact result for the bare Argonne V18 interaction and the 
experimental binding energy are indicated by horizontal lines. 



interactions play a somewhat different role. In the case of the UCOM(var.) interactions 7,^ only 
affects the range of the tensor correlator in the deuteron channel. In the case of UCOM(SRG) and 
SRG interactions a affects also the central correlations in all channels. Beyond a certain point, 
stronger central correlations actually result in less binding from the central part of the interaction 
and the interaction becomes less attractive. 

4.2. Tjon-line and the role of three-body interactions 

As has already been observed by Tjon for local interactions, a correlation exists between the 



and "^He binding energies [17 311 . When the binding energy of "^He is plotted against the binding 
energy of the results for different interactions fall essentially onto a single line, the so-called 
Tjon-line. This has been confirmed also for modern interactions. Typically the binding energies 
for bare two-body forces are too small compared to experiment. This can be corrected by adding 
an appropriate three-body force. In Fig.[2l]results from Refs. [10, s]] for realistic two-body interac- 



tions and combinations of two- plus three-body forces are shown. In addition we show the results 
with UCOM and SRG interactions for which the binding energies vary as a function of tensor 
correlation range or flow parameter a and continue more or less on the Tjon-line obtained from 
bare realistic two-body forces. The results obtained with the UCOM and SRG interactions are 
very similar and the trajectories of all interactions pass very close (within 200 keV) to the exper- 
imental binding energies of and '^He. The results closest to the experimental binding energies 
are obtained for the parameters: 

UCOM(var.): /^ = 0.09fm^ 

UCOM(SRG): a = 0.04fm^ (125) 
SRG: ar = 0.03fm'* 

This is a remarkable result. By choosing a particular tensor correlation range or particular 
flow parameters the two-body UCOM or SRG interactions not only reproduce the experimental 
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Figure 21: Binding energy of "^He versus binding energy of calculated with UCOM(var), UCOM(SRG) and SRG 
interactions. Also included are results with bare two-body forces (black diamonds) and combinations of two- and 
three-body forces (grey diamonds) taken from lilQil and 



nucleon-nucleon scattering data but also give the correct binding energies in the three- and four- 
body systems. In this calculation we neither evaluate three- and four-body contributions from the 
correlated two-body interaction nor do we include genuine three- and four-body forces. At this 
point all three- and four-body contributions have to cancel at least on the level of the expectation 



value. As explained by Polyzou and Glockle Ill4ll a unitary transformation exists between different 



combinations of on-shell equivalent two- and three-body interactions. UCOM and SRG interac- 
tions provide a particular realization of such a transformation. The three-body contributions of 



SRG interactions have been studied explicitly in m% by evolving the three-body matrix elements 
in the harmonic oscillator basis. 

We will use UCOM and SRG interactions with tensor correlation range and flow parameters 
optimized for three- and four-body systems as described above for NCSM calculations in p-shell 
nuclei in Sec. 14.31 and for exploratory studies using Hartree-Fock and many-body perturbation 
theory for doubly magic nuclei up to ^"^Pb in Sec. [5l 

4.3. Properties of A = 6,1 nuclei 

In this section we study properties of the light ju-shell nuclei ^He, '^Li and ^Li using the NCSM 
code Antoine ll74l ItsIi . Matrix elements of UCOM and SRG interactions are provided as jj- 
coupled matrix elements generated from the relative matrix elements using the Talmi-Moshinsky 
transformation. Again we do not employ the Lee-Suzuki transformation. 

Like for the three- and four-body systems the results are studied as a function of oscillator 
frequency hQ. and model space size A^max- We do not reach full convergence in the accessible model 
spaces for these nuclei, even with the UCOM or SRG interactions. For ground-state energies we 
can use extrapolations to estimate the converged binding energy. An error estimate is provided 
by comparing results from different oscillator parameters. For observables like the radii and the 



55 





UCOM(var.) 


UCOM(SRG) 


SRG 


Experiment 


"He 


27.9(4) 


28.4(3) 


28.8(4) 


29.269 




30.9(4) 


31.6(3) 


32.0(4) 


31.995 




37.4(6) 


38.7(4) 


39.6(5) 


39.245 



Table 3: Experimental and extrapolated calculated binding energies (in MeV) obtained with UCOM(var.), 
UCOM(SRG) and SRG interactions. Error estimates are obtained by comparing extrapolated energies for different 
oscillator frequencies. 



quadrupole moment of ^Li the results show a strong dependence on the model space size and 
extrapolated results are not reliable. Spectra on the other hand appear to be much better converged 
and can be compared with experiment. 

The convergence problems are to a large extend not caused by properties of the interactions 
but by the cluster or halo nature of these nuclei. To properly describe the asymptotics of wave 
functions with a neutron halo like in "He, or an underlying cluster structures in case of the lithium 
isotopes, large model spaces are needed in the oscillator basis. 

Ground-state energies 

We calculate the ground-state energies in model spaces up to A^^ax = 14 for ^He and ^Li and 
up to A^max = 12 for ^Li for oscillator frequences HQ. = 12, 16, . . . , 28 MeV. For all interactions the 
lowest energy in the OhQ space is obtained for an oscillator frequency of 16 MeV. In the largest 
model spaces lowest energies are found for oscillator constants between 24 MeV and 28 MeV for 
the UCOM(var.) interaction and between 20 MeV and 24 MeV for UCOM(SRG) and SRG inter- 
actions. To estimate the converged ground-state energy an exponential extrapolation is performed 
using the results obtained in the four largest model spaces. An estimate of the error can be obtained 
by comparing the extrapolated energies for different oscillator frequencies. In Figs. I22l23l24l the 
calculated energies are shown as a function of the model space size A^max together with the fitted 
exponentials. Both UCOM(SRG) and SRG interactions provide extrapolated binding energies that 
are close to the experimental values, whereas the binding energies with the UCOM(var.) interac- 
tion are somewhat underestimated for all nuclei. The results are summarized in Table [3l Note 
that we use the interaction parameters (11251) as obtained from the Tjon-line analysis — no further 
adjustments are made here or in the following. 

Radii, magnetic dipole moments and quadrupole moments 

Electromagnetic properties provide important tests of the wave function beyond the simple 
binding energy. We calculate the point proton radii for all nuclei and the magnetic dipole moment 
as well as the electric quadrupole moment for ^Li and ^Li. The results are very similar for the 
different UCOM and SRG interactions and we only show results for the UCOM(SRG) interaction. 
The results are also very similar to results from NCSM calculations using the CD Bonn and INOY 



interactions 11761. IVVH. 

The point radii shown in Fig. [25] are calculated for the oscillator frequencies hD. = 16, 20, and 
24 MeV. There still is a strong dependence of the calculated radii on the model space size. We 
do not believe that the extrapolation of the radii is reliable, but the results indicate that the radii 
are slightly too small when compared to experiment. It is difficult to draw conclusions regarding 
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Figure 22: Energy of the ^He 0^ state as a function of model space size for different oscillator frequencies obtained 
with UCOM(var.), UCOM(SRG) and SRG interactions. Exponential extrapolations are fitted to the results from the 
four largest model spaces. 
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Figure 23: Energy of the ^Li 1^ state as a function of model space size for different oscillator frequencies obtained 
with UCOM(var.), UCOM(SRG) and SRG interactions. Exponential extrapolations are fitted to the results from the 
four largest model spaces. 
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Figure 24: Energy of the ^Li 3 /2" state as a function of model space size for different oscillator frequencies obtained 
with UCOM(var.), UCOM(SRG) and SRG interactions. Exponential extrapolations are fitted to the results from the 
four largest model spaces. 
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Figure 25: Point proton radii of ^He, ^Li and ^Li as a function of model space size for different oscillator frequencies 
using the UCOM(SRG) interaction. Experimental values from 117 811 and 17911 . 
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Figure 26: Magnetic dipole moment of *Li and 'Li as a function of model space size calculated with UCOM(SRG). 
Experimental values from Isoll . 



the saturation properties of the interaction as the radii for these nuclei depend strongly on the 
asymptotic behavior of the wave function due to their halo or cluster nature. 

The magnetic dipole moments of ^Li and 'Li shown in Fig. [26] agree reasonably well with the 
experiment. The quadrupole moment of ^Li is correctly predicted to be very small and negative, 
which confirms the "^He plus deuteron picture. For the quadrupole moment of 'Li we find a sim- 
ilar behavior as for the radii. There is still a strong dependence on the model space size. The 
extrapolated results underestimate the experimental value. 

Spectra 

Further information about the interaction can be obtained from the spectra of excited states. 
The spectra are calculated in llhD. and lOhQ model spaces for the A = 6 and A = 7 nuclei 
respectively. We show here the results for the oscillator parameter 16 MeV which corresponds to 
the ground state minimum in the OhD. model space and which shows the fastest convergence for 
the spectra. For other oscillator constants the spectra vary more rapidly when enlarging the model 
space but the converged results depend only very weakly on the oscillator parameter. 

The results are very similar for all three interactions. In Fig. [28] the results are summarized. 
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Figure 27: Quadrupole moment of ^Li and ^Li ground state as a function of model space size for UCOM(SRG). 
Experimental values from Isoll . 



For the UCOM(var.) interaction the results are shown starting from the OhQ. model space. For the 
UCOM(SRG) and SRG interactions we only show the results from the largest model spaces. 

In ^He the energy of the 2^ state is well converged. This is not true for the second 2^ and the 1"^ 
state. These states are well above the two-neutron separation energy and there is no experimental 
confirmation for the existence of these states. 

For ^Li the T = states are well converged in contrast to the T = \ states (22 and 0^). We 
also find the excitation energy of the 3^ state to be too high. This indicates that the eff"ective spin- 
orbit force in the UCOM and SRG interactions is too weak. A similar observation can be made in 
the spectrum of 'Li. Here the splittings between the 3 /2" and 1 /l" states as well as between the 
7 /2~ and 5/2" states are too small compared to the experimental values, as well. In that respect 
the UCOM and SRG interactions perform similarly to other two-body interactions. It has been 
observed in GFMC calculations Ul\\ that three-body forces contribute significantly to the eff"ective 
spin-orbit strength. NCSM calculations with chiral two- and three-body forces [18 ill also show a 
strong dependence of these splittings on the parameters of the three-body force. 
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Figure 28: Spectra of low lying states in ''He, <'Li and ^Li calculated with UCOM(var), UCOM(SRG) and SRG 
interactions in comparison with experiment Isoll . 
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5. Hartree-Fock and beyond 



For nuclei beyond the p-shell, many-body calculations in methods like the full NCSM are 
not feasible anymore, because the dimension of the many-body basis at a given N^ax^^ truncation 
level grows factorially with the particle number. One can extend the domain of NCSM calculations 
to larger particle numbers and model space sizes by using importance truncation methods, as 



discussed in Refs. i82Ll83l . l84|] . For the description of ground states of closed-shell and neighboring 



nuclei, coupled cluster methods have been employed quite successfully l.85i.i86ui87i.i88i . i89i.i9Q . i9in . 
These methods are able to cover the majority of correlations in the nuclear many-body system, 
however, they are computationally demanding as well. 

At the opposite end of the scale regarding the computational cost and the ability the describe 
correlations is the Hartree-Fock approach ll92ll . A simple Hartree-Fock calculation — based on 
a single Slater determinant for the description of the ground state — can be done easily for any 
isotope throughout the nuclear chart. Obviously, the Hartree-Fock approach does not allow for 
the description of any correlations and, therefore, cannot provide a quantitative approximation 
for nuclear observables when using realistic nuclear interactions. It does, however, provide a 
variational upper bound for the exact ground-state energy and as such can be used to assess the 
qualitative systematics, e.g., of the binding energies as function of mass number, throughout the 
whole nuclear chart. Furthermore, the Hartree-Fock solution can serve as a starting point for 
improved approximations that take the missing correlations into account. In the simplest case, 
low-order many-body perturbation theory can be used to estimate the effect of correlations on 
the energy or other observables. Alternative methods for the inclusion of correlations beyond 
Hartree-Fock, such as ring-, ladder-, or Pade-resummed perturbation theory 11931 19411 as well as 
Brueckner-Hartree-Fock schemes and Green's function methods [1951 . 19611 are also feasible. The 
Hartree-Fock or the corresponding Hartree-Fock-Bogoliubov 11971 . 19811 solutions also form the basis 
for the descr iption of collective excitations in the Random Phase Approximation at different orders 

In the following we will use the Hartree-Fock scheme as well as second-order many-body 
perturbation theory to study the systematics of binding energies and charge radii resulting from 
the UCOM- and SRG-transformed interactions. The aim is not to provide a precise prediction of 
the ground state energy for heavy nuclei, but to assess the systematic behavior of the transformed 
two-body interactions with increasing mass number. 



5.1. Hartree-Fock with UCOM- and SRG-transformed interactions 

In the simplest formulation of a Hartree-Fock (HF) scheme, the many-body state is approxi- 
mated by a single Slater determinant 

|HF) = |Om) = ( 10,,) ® 10,,) ® • • • ® 10,,)), (126) 

where Jl is the antisymmetrization operator acting on an A-body product state. The single-particle 
states |0y) are used as variational degrees of freedom in a minimization of the expectation value 
of the many-body Hamiltonian. The formal variational solution of the many-body problem using 
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the trial state (11261) leads to the w ell kn own HF equations for the single-particle states, which have 



to be solved self-consistently i92Lll0in . 

The Hamiltonian itself is the same as it was used in Sec. |4]for the NCSM calculations. It 
consists of the intrinsic kinetic energy Tint = T - T^m and the transformed two-body interaction 
Vnn including Coulomb and charge-dependent terms 

Hint = T - Tcm + Vnn = Tint + Vnn > (127) 

Unlike the NCSM, the use of this translational invariant Hamiltonian does not guarantee that the 
HF ground state is free of spurious center-of-mass contaminations. The Slater determinant form of 
the many-body state does not allow a separation of intrinsic and center-of-mass motion for general 
single-particle states. However, for the purpose of the present discussion, the effect of center-of- 
mass contaminations on the ground-state energy is irrelevant. A more stringent but computation- 
ally expen sive treatm ent of the center-of-mass problem would require an explicit center-of-mass 



projection Ill0lll03h 



We formulate the HF scheme in a basis representation using harmonic oscillator single-particle 
states. Thus the matrix elements entering HF equations are the same as in the NCSM calculations 
of Sec. |4j The HF single particle states |0v) are written as 

10,) = \aljmm,) = J] '"""'^ \nljmm,) , ^28) 

n 

where \nl jmnit) denotes a harmonic oscillator single-particle state with radial quantum number n, 
orbital angular momentum /, total angular momentum j with projection m, and isospin projection 
quantum number m,. Assuming spherical symmetry, only oscillator states with the same quantum 
numbers /, j, and m can contribute in the expansion. In the following, we will restrict ourselves 
to constrained or closed-shell calculations, where c^"^^"™'^ = Cji"'^'"'^ is independent of m. The 
details of the resulting HF equations for the expansion coefficients C'^"'-'"''^ and of their solution are 



discussed in Ref. lllOlh 



Within the HF approximation, we consider the ground-state energies and the charge radii for 
a sequence of nuclei with closed j-shells from "^He to ^°^Pb using the different UCOM- and SRG- 
transformed interactions adopted in Sec. IH For these calculations the harmonic oscillator single- 
particle basis includes 15 major shells, i.e., the sum in (11281) is limited to 2n -i- / < 14, which is 
sufficient to guarantee full convergence of the HF states for all nuclei under consideration. We 
use a sequence of oscillator lengths ^hg from 1.3 fm to 2.4 fm for the underlying oscillator basis. 
In accord with the variational principle, we adopt the value of the oscillator length that yields the 
minimal ground-state energy, though, for the basis sizes used here, the energies and radii obtained 
on the HF level are largely independent of the oscillator length. 

A summary of the HF ground- state energies and charge radii obtained with the UCOM and 
SRG-transformed potentials is given in Fig. [29j Here we adopt the UCOM-correlators obtained 
variationally and through the SRG-mapping with the range or flow parameters determined by 
fitting to the experimental "^He ground-state energy in converged NCSM calculations. Similarly, 
the flow parameter for the SRG-evolved interaction is determined through the "^He binding energy 
(cf. Sec. 14.21) . The HF approximation, therefore, yields similar results for the ground-state energy 
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Figure 29: Ground-state energies per nucleon (a) and charge radii (b) for a sequence of nuclei with closed j-shells 
obtained at the Hartree-Fock level. The data sets correspond to UCOM(var) interaction (•), the UCOM(SRG) 
interaction (■), and the SRG interaction (♦). The range or flow parameter of the different transformations is fixed 
such that the experimental "*He ground-state energy is reproduced in NCSM calculations (cf. Sec. 14.21 ). The black 
bars indicate experimental values ll04ll . 



of "^He with all three interactions — UCOM(var.), UCOM(SRG), and SRG. However, the HF energy 
is above the experimental and the converged NCSM ground-state energy, which is expected — the 
single determinant describing the HF ground state corresponds to a Oh£l NCSM eigenstate and 
cannot describe any of the correlations that the NCSM model space will capture with increasing 
model space size Nmax^^- 

With increasing mass number the behavior of the UCOM and the SRG-transformed interac- 
tions differs dramatically. The UCOM interactions, both UCOM(var.) and UCOM(SRG), lead to 
HF ground- state energies per nucleon that are of the order of 4 or 5 MeV throughout the whole 
mass range. Especially for the UCOM(var.) interaction there is a practically constant offset of 
about 4 MeV per nucleon between the HF energies and experiment. For the case of "^He we know 
from the NCSM calculations in Sec. |4]that the missing binding energy resides in correlation en- 
ergy, because the same effective interaction reproduced the experimental values in large model 
spaces. We expect that the inclusion of correlations beyond HF will also lower the ground-state 
energies of the heavier nuclei and bring them closer to experiment. 

The SRG-transformed interaction exhibits a very different trend: The HF binding energy per 
nucleon increases rapidly with increasing mass number. Already for intermediate masses, the HF 
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energy drops below the experimental ground-state energy. Given that the HF energy gives an upper 
bound for the exact energy eigenvalue of the Hamiltonian, this discrepancy cannot be remedied 
through the inclusion of beyond-HF correlations, but hints at the induced many-body forces (see 
discussion in Sees. [Tl and [231 that are left out in this calculation. In the SRG they have an over-all 
repulsive effect. 

This intrinsic difference between UCOM and SRG interactions is also reflected in the charge 
radii depicted in Fig. [29l b). The charge radii obtained with the UCOM interaction show a sys- 
tematic deviation form the experimental trend. The predicted radii are too small and the difference 
to experiment increases linearly with increasing mass number reaching a deviation of about 1 fm 
for ^°^Pb. For the SRG interaction, the deviation is even more pronounced, the radius of ^o^pj, 
underestimated by about 2 fm. 

There are two main differences between the UCOM-transformed and the SRG-transformed 
two-body interactions which cause the different behavior when going to heavier nuclei: (?) The 
UCOM interactions use a unitary transformation that is optimized to account for short-range cor- 
relations in the lowest partial wave of each spin-isospin channel only. Thus, higher partial waves, 
whose impact grows with increasing mass number, are not pre-diagonalized in an optimal way 
by the UCOM transformation. However, because the wave functions for larger L are suppressed 
at short distances by the centrifugal barrier, this is a minor effect. The SRG transformation, in 
contrast, handles each partial wave separately and thus leads to an optimal pre-diagonalization 
for all. (zz) Even for the lowest partial waves the UCOM transformation does not provide the 
same perfect pre-diagonalization in the high-momentum sector as the SRG transformation, as dis- 
cussed in Sec. 13.71 The residual off-diagonal high-momentum matrix elements together with the 
less pronounced pre-diagonalization of the higher partial waves stabilize the UCOM interactions 
against the overbinding observed in the SRG calculations. On the other hand this leads to slower 
convergence of the UCOM interactions as compared to the SRG. 

One can view the difference between SRG and UCOM from yet another perspective. Both 
approaches use a unitary transformation, which preserves the eigenvalues of the Hamiltonian in 
many -body space, provided the transformation is done without any additional truncations. How- 
ever, here we use the cluster expansion and truncate at two-body level, i.e., we discard all the 
induced many-body forces. Thus, if we observe a systematic difference in an exact many-body 
calculation using two-body part of the transformed interactions only, we can conclude that the 
omitted many-body forces behave systematically different in SRG and UCOM, because the trans- 
formed Hamiltonians in their complete form in A-body space have to yield the same result. In 
this view, the simplistic HF calculations already show that the induced three-body and many-body 
interactions in the SRG framework must have a much larger net effect on binding energies and 
wave functions — and thus on radii — than in the UCOM framework. 



5.2. Low-order many-body perturbation theory 

The simplest way to estimate the effect of correlations beyond HF is low-order many-body 
perturbation theory. Many-body perturbation theory (MBPT) starting from the HF solution is a 
stand ard te chnique in many fie lds of quanturn many-body ph ysics, ranging from quantum chem- 



istry 111 0511 to nuclear physics lll06l 1107L llOSL 1109L llOlL l94ll . It is straightforward to apply and 
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^He 40ca 48Ni eo^i 885^ 100s„ 132s„ 208pb 

I60 34si 48ca 56Ni 78Ni 90^^ lU^^ U6q^ 



Figure 30: Ground-state energies for a sequence of nuclei with closed j-shells obtained at the HF level (open symbols) 
and in second-order MBPT (filled symbols) using the UCOM(var.) interaction. The different curves correspond to 
different model spaces built from 11 (•), 13 (■), and 15 (♦) major oscillator shells. The oscillator length for each 
nucleus is chosen according to the HF root-mean-square radius (a) or according to the experimental charge radius (b). 
The black bars indicate the experimental binding energies 110411 . 



computationally simple, but has inherent limitations. It is well known that the convergence of sue 



cessive orders of perturbation theory is not guaranteed, on the contrary 111 101 . 1941] . Nevertheless, 
low-order MBPT provides at least a qualitative measure for the effect of correlations beyond HF. 

We will restrict ourselves to second-order calculations for an order-of-magnitude estimate of 
the correlation energy, i.e., the change in the ground-state energy resulting from beyond-HF cor- 
relations. The second order contribution involves anti symmetrized two-body matrix elements of 
the intrinsic Hamiltonian Hint containing two HF single-particle states below the Fermi energy 
(denoted by v, v') and two HF single-particle states above the Fermi energy (denoted by iu,fi'y. 



K0v0v'|Hint|0^0^')P 



V.V' fl.^t' I- , fA )^ 



(129) 



Note that the full two-body part of the many -body Hamiltonian enters, which includes the intrinsic 
kinetic energy in our case. 

The starting point for the evaluation of the correlation energy via perturbation theory is the 
HF solution yielding a finite set of single-particle states |0v) and the corresponding single-particle 
energies for the respective nucleus. As discussed before, we use the oscillator basis including a 
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^He ^C^ 48Ni eo^i 885^ ioos„ i32s„ 208pb 

I60 34si 48ca 56Ni 78Ni 90^^ lUgj^ 146^^ 

Figure 3 1 : Ground-state energies for a sequence of nuclei with closed j-shells obtained at the HF level (open symbols) 
and in second-order MBPT (filled symbols) using the UCOM(SRG) interaction. The different curves correspond to 
different model spaces built from 11 (•), 13 (■), and 15 (♦) major oscillator shells. The oscillator length for each 
nucleus is chosen according to the HF root-mean-square radius (a) or according to the experimental charge radius (b). 
The black bars indicate the experimental binding energies 110411 . 



certain number of major shells for a specific oscillator length. Because perturbation theory is not a 
variational approach, we cannot use variational arguments to find an optimal oscillator length, but 
have to resort to other prescriptions. We adopt two dilFerent schemes for choosing the oscillator 
length aHo for each nucleus: Using either the root-mean- square radius predicted by the HF solution 
or the experimental charge radius we optimize the oscillator length such that a naive shell-model 
Slater determinant built from harmonic oscillator single-particle states approximately reproduces 
the respective radius. Because the HF solutions for the UCOM interaction underestimate the 
charge radii, the oscillator lengths obtained from the HF radius are smaller than the ones obtained 
from the experimental radius. 

The ground-state energies obtained by including the second-order MBPT contribution on top 
of the HF energy are presented in Figs. [30] and [3T] for the UCOM-transformed interactions with the 
UCOM(var.) and the UCOM(SRG) correlators, respectively. The upper and lower panels in each 
figure are obtained using the HF radii and the experimental radii, respectively, to fix the oscillator 
length. For each case we show a sequence of calculations in model spaces consisting of 11, 13, 
and 15 major oscillator shells in order to assess the convergence behavior. The HF energies are 
fully converged and independent of the model space size. The second-order MBPT estimate of 
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the correlation energy (11291) shows signatures of convergence only for light isotopes, for heavier 
isotopes there still is a significant change of typically 0.5 MeV per nucleon when going from 13 to 
15 shells. Uncertainties of a similar order of magnitude result from the dependence of the second- 
order energy on the oscillator length, with decreasing auo the second-order energy contribution 
l^^^^l is increasing. Finally, one should keep in mind that second-orderJMBPT provides only a 
crude approximation for the correlation energy. As shown in Ref. [i94ll the deviations of the 
second-order estimate to the exact eigenvalue in the same model space can be sizable. 

Despite the uncertainties associated with the second-order MBPT calculation regarding con- 
vergence, choice of the oscillator length, and quality of low-order MBPT as such, the results in 
Figs. |30] and |3T] prove that the correlations beyond HF cause an essentially constant shift of the 
ground state energy per nucleon across the whole mass range. For the UCOM interactions, this 
brings the ground state energies into the same regime as the experimental binding energies. The 
over-all systematics of binding energies obtained with the UCOM interactions is in agreement with 
experiment, already at the level of a pure two-body interaction. The role of three-body interactions 
is reduced to providing corrections, e.g., regarding the charge radii, on top of an already reason- 
able trend. For the standard SRG-transformed interactions the situation is different: The two-body 
component alone cannot provide the correct systematics and three-body interaction have to have a 
strong impact on the binding-energy systematics already. 
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6. Conclusion 



The Unitary Correlation Operator Method (UCOM) provides a universal tool to account for 
short-range interaction-induced correlations in the nuclear many-body problem. The correlations 
can be either imprinted into many-body states that otherwise would not contain these correlations, 
like Slater determinants, or they can be absorbed in effective operators that are defined through an 
explicit similarity transformation. The effective interactions obtained in this way are well suited 
for low-momentum Hilbert spaces, because, unlike bare interactions, they do not scatter strongly 
to high momenta. UCOM is very transparent and intuitive as it explicitly introduces correlation 
functions for the description of short-range central and tensor correlations. These correlation 
functions play the role of variational degrees of freedom for the many-body states. We propose 
two methods to find optimal correlation functions for a given bare Hamiltonian: One is energy 
minimization in two-body space; the other employs the close relationship to the SRG approach, 
which also aims at separating low- and high-momentum scales. 

UCOM has been developed for both, matrix representation and operator representation. The 
matrix representation in the harmonic oscillator basis is used for NCSM calculation^. We show 
that binding energies and spectra converge much more rapidly with increasing size of the Hilbert 
space when using effective UCOM-transformed interactions rather than the bare interaction. It 
is even possible to do Hartree-Fock calculations and obtain bound nuclei throughout the nuclear 
chart. We show that this ab initio Hartree-Fock method can only account for about half the binding 
energy, the other half is correlation energy that cannot be obtained by a single Slater determinant. 

Especially the operator representation of UCOM explains in a very transparent way, why strong 
short-range correlations and long mean free path or the mean-field shell model are not contradict- 
ing each other. There is a separation of length scales, well developed for the central correlations 
and less well for the tensor correlations, which allows to renormalize the bare Hamiltonian to 
an effective one appropriate for low-momentum Hilbert spaces. One might be tempted to be- 
lieve that short-range correlations are not real, as soft phase-shift equivalent interactions with a 
modified off-shell behavior describe the asymptotic properties of the two-nucleon system and are 
successfully used for the description of many-body systems. However, short-range correlations 
are revealed by the stiffness of nuclear matter against compression or the high-momen tum tails o f 
nucleon momentum distributions that are presently gaining renewed interest [1481 pOl IlllL 111211 . 
Around the Fermi edge the momentum distributions differ from the mean-field ones also due to 
long-range correlations that are not accounted for by UCOM. These correlations have to be treated 
by configuration mixing. 

As the short-range correlations are, to a large extend, state independent and can be treated 
by a unitary transformation, one can work in the independent-particle basis using a transformed 
Hamiltonian. Also electric and magnetic observables are seeing mainly the low-momentum nature 
of the states and change little by the UCOM transformation. An exception are Gamow-Teller 
transitions, which are sensitive to tensor correlations. It is well known that a quenching of typically 
0.8 occurs for Hilbert spaces that do not contain these correlations. 



-The UCOM two-body matrix elements in the harmonic oscillator basis (relative LS - or yy-coupling for various 
frequencies and model-space sizes) suitable for no-core calculations are available from the authors upon request. 
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It turns out that effects from UCOM-induced three-body interactions cancel to a certain extend 
those from the original three-body force. The partial cancellation effect is not yet understood. 
Heavier nuclei and the saturation properties of nuclear matter indicate that three-body forces can- 
not be substituted completely. The short-range repulsion of the nuclear interaction is essential for 
describing the correct saturation properties. In particular when increasing the density above twice 
nuclear saturation density, the short-range repulsive correlations are expected to become so strong 
that they cannot be treated anymore by a two-body approximation for the effective interaction. 
From the above arguments it is clear that any ab initio treatment of nuclear matter at higher densi- 
ties based on Slater determinants of single-particle plane waves demands a sophisticated effective 
interaction with many -body forces. 

The Unitary Correlation Operator Method has been formulated in momentum representation, 
in harmonic oscillator basis and in operator representation. The latter can be used in any represen- 
tation for example in Fermionic Molecular Dynamics (FMD), which provides many-body Hilbert 
spaces especially suited for cluster structures in nuclei (not discussed in this paper). UCOM is a 
very versatile approach that provides effective interactions as well as the corresponding effective 
operators. 
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